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Abstract 

We evaluate all two-body decay modes of the heavy scalar top quark in the Minimal Su- 
persymmetric Standard Model with complex parameters (cMSSM) and no generation 
mixing. The evaluation is based on a full one-loop calculation of all decay channels, 
also including hard QED and QCD radiation. The renormalization of the complex 
parameters is described in detail. The dependence of the heavy scalar top quark decay 
on the relevant cMSSM parameters is analyzed numerically, including also the decay 
to Higgs bosons and another scalar quark or to a top quark and the lightest neutralino. 
We find sizable contributions to many partial decay widths and branching ratios. They 
are roughly of 0(10%) of the tree-level results, but can go up to 30% or higher. These 
contributions are important for the correct interpretation of scalar top quark decays 
at the LHC and, if kinematically allowed, at the ILC. The evaluation of the branch- 
ing ratios of the heavy scalar top quark will be implemented into the Fortran code 
FeynHiggs. 
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1 Introduction 



One of the most important tasks at the LHC is to search for physics effects beyond the 
Standard Model (SM), where the Minimal Supersymmetric Standard Model (MSSM) [1] 
is one of the leading candidates. Supersymmetry (SUSY) predicts two scalar partners for 
all SM fermions as well as fermionic partners to all SM bosons. Another important task 
is investigating the mechanism of electroweak symmetry breaking. The most frequently 
investigated models are the Higgs mechanism within the SM and within the MSSM. Contrary 
to the case of the SM, in the MSSM two Higgs doublets are required. This results in five 
physical Higgs bosons instead of the single Higgs boson in the SM; three neutral Higgs 
bosons, hn [n = 1,2,3), and two charged Higgs bosons. 

If SUSY is realized in nature and the scalar quarks and/or the gluino are in the kine- 
matic reach of the LHC, it is expected that these strongly interacting particles are copiously 
produced. This includes the production of scalar top quark pairs or the production of two 
gluinos with the subsequent (possible) decay to a scalar top quark and a top quark. An 
interesting production channel of Higgs bosons at the LHC is the decay of the heavy scalar 
top quark to the lighter scalar top (scalar bottom) quark and a neutral (charged) Higgs 
boson, see, for instance, Refs. [2l[3] and references therein. At the ILC (or any other future 
e~^e~ collider such as CLIC) a precision determination of the properties of the observed par- 
ticles is expected [HIS]. (For combined LHC/ILC analyses and further prospects see Ref. [6].) 
Thus, if kinematically accessible, Higgs production via scalar top quark decays could offer 
important information about the stop and Higgs sector of the MSSM. 

In order to yield a sufficient accuracy, one-loop corrections to the various scalar top quark 
decay modes have to be considered. We take into account all two-body decay modes of the 
heavy scalar top quark, ^2, in the MSSM with complex parameters (cMSSM), but we neglect 
flavor violation effects and resulting decay channels that rather play a role for the decay of 
the light scalar top quark, ti, in special regions of the MSSM parameter space [7]. More 
specifically, we calculate the full one-loop corrections to the partial decay width^ 



rik^hhn) (n = 1,2,3), (1) 

r(t2 ^ hZ) , (2) 

rih^txl) (fc = 1,2,3,4), (3) 

rih ^ t~g) , (4) 

r{i2^kH+) {t = l,2), (5) 

r{i2^kW+) (z = l,2), (6) 

ni2^bxt) (J = l,2), (7) 



where Xk denotes the neutralinos, g the gluino, xf the charginos, t and b the top and bottom 
quark and Z and the SM gauge bosons. The total decay width is defined as the sum of 
the partial decay widths ([I]) to ([7]), where for a given parameter point several channels may 
be kinematically forbidden. 

^ It should be noted that the purely loop induced decay channels ^2 hi/ 9 have been neglected because 
they yield exactly zero; see Sec. [3] for further details. 
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As explained above, we are especially interested in the branching ratios (BR) of the decays 
of the ^2 to a Higgs boson and another squark, Eqs. ([T]) and (jSj), as part of an evaluation 
of a Higgs production cross section. This can be an interesting production channel at the 
LHC (see, for instance, Ref. |8], where pp — )■ i\iih is analyzed, or Ref. [9], where searches for 
CP-odd Higgs bosons in top squark decays are discussed). However, in order to reach a high 
accuracy, all two-body decay channels should be evaluated at one-loop. On the other hand, 
because we are interested in two-body modes (involving Higgs bosons), it is not necessary 
to investigate three- or four-body decay modes as these only play a significant role once the 
two-body decay modes are kinematically forbidden, and thus the relevant BR are zero. 

We also concentrate on the decays of ^2 and do not investigate 4 decays. In the presence 
of complex phases this would lead to somewhat different results. However, such an analysis 
of CP-violating effects is beyond the scope of this paper. 

Higher-order contributions to scalar fermion decays have been evaluated in various analy- 
ses over the last decade. However, they were in most cases restricted to one specific channel. 
In many cases only parts of a one-loop calculation has been performed, and no higher-order 
corrections in the cMSSM are available so far. More specifically, the available literature 
comprises the following. First, 0{as) corrections to partial decay widths of various squark 
decay channels in the MSSM with real parameters (rMSSM) were derived: to the decay of 
a squark to a quark and a chargino or neutralino in Ref. [10], to the decay of a squark to a 
quark and a gluino in Refs. [II1[I2], to the decay of a squark to a squark and a SM gauge 
boson in Ref. [13] , and to the neutral Higgs boson radiation of a scalar top or bottom quark 
in Ref. [13]. (Those 0{as) corrections have been implemented into the code SDECAY p5].) 
A tree- level analysis on several t and b decay modes was presented in Ref. [TB]. In a second 
step, Yukawa corrections to the partial decay widths of a scalar quark were evaluated in 
Refs. [T7l[18] . Finally, full one- loop contributions were derived, for the decay of a squark to 
a quark and a chargino or neutralino in Ref. [19], and for the decay of a scalar fermion to a 
scalar fermion and a gauge boson in Ref. [20]. One-loop corrections to scalar quark decays 
in the rMSSM, derived in a pure DR scheme (see below) have been made available in the 
program package SFOLD [21]. Also the partial decay width of a CP-even and a CP-odd Higgs 
boson to scalar quarks at the one-loop level is available; see Refs. [22] and [23], respectively. 
A more recent evaluation can be found in Ref. [21]. Tree- level analyses for the decay of a 
t or a 6 in the cMSSM have been published in Refs. [251126] . and a one-loop calculation of 
a scalar top quark decaying to a bottom quark and a chargino in the cMSSM is presented 
in Ref. [27], where an LHC specific analysis can be found in Ref. [2H]. Finally, results in an 
effective Lagrangian approach can be found in Ref. [29]. 

Several methods have been discussed in the literature to extract the complex parame- 
ters of the model from experimental measurements. Branching ratios at a linear collider 
were analyzed at the tree- level in Ref. [2^j. Triple products of decaying scalar top or bot- 
tom quarks have been examined in Ref. [SU] (without specifying the production modes) and 
in Refs. [3TI433] at the LHC. Rate asymmetries for decaying top squarks are analyzed in 
Ref. [27], again without specifying the production modes, and especially for the LHC in 
Ref. [28]. Depending on the realized cMSSM parameter space and on some further assump- 
tions on the LHC performance, it seems to be possible to obtain limits on, e.g., the phases 
of Ml, At and Af, at the LHC. No corresponding analysis, to our knowledge, of the phase of 
M3 has been performed so far. 
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In this paper we present for the first time a full one-loop calculation for all two-body 
decay channels of the heavier scalar top in the cMSSM (with no generation mixing), taking 
into account soft and hard QED and QCD radiation. In Sec. [2]we review the renormalization 
of all relevant sectors of the cMSSM. Details about the calculation can be found in Sec. [3l 
and the numerical results for all decay channels are presented in Sec. HI The conclusions can 
be found in Sec. [51 The results will be implemented into the Fortran code FeynHiggs [3H437] . 

2 The complex MSSM and its renormalization 

All the channels ([1]) - ([7j) are calculated at the one-loop level, including hard QED and QCD 
radiation. This requires the simultaneous renormalization of several sectors of the cMSSM, 
including the colored sector with top and bottom quarks and their scalar partners as well as 
the gluon and the gluino, the Higgs and gauge boson sector with all the Higgs bosons as well 
as the Z and the W boson and the chargino/neutralino sector. In the following subsections 
we briefly review these sectors and their renormalization. To our knowledge, it is the first 
time that such a complete renormalization of the cMSSM has been performed. 

2.1 The colored sector of the cMSSM 

The colored sector of the cMSSM can be divided into a quark/squark part, a gluino part 
and a gluon part. The quark/squark part contains the soft SUSY-breaking mass parameters 
Mq^, Mg^, the trilinear couplings Ag, the quark masses rnq as well as the quark and the 
squark field^ q and g, while the gluino part comprises the soft SUSY-breaking gaugino mass 
parameter M3 and the gluino field g. From the gluon part, only the renormalization of the 
strong coupling constant is needed for our calculation. 

2.1.1 The top and bottom quark/squark sector 

The part of the Fourier transformed Lagrangian that is bilinear in the quark and the squark 
fields with q = {t, b} and q = {i, b} can be written as 

= {ql (P'li - M,-) (^J^^ + - m,)u.q + - m,)u+q , (8) 

where u± = ± 75) are the right- and left-handed projectors, respectively, iriq with 
q = {t, b} is the {top, bottom} quark mass and the stop and sbottom mass matrices, M^- 
and Mj^, are given by 

' V ^,X, Ml+ml + Mlc2^Q,sl 

with 

Xq = Aq—fi*hi, K = {cot /3, tan/3} for q = {t,b}, 02/3 = cos 2/3. (10) 

^ It should be noted that for the renormahzation of the quark/squark sector we focus on the thhd 
generation - which is the relevant part for our calculation - but, in principle, it can be generalized to the 
other generations. 
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The soft SUSY-breaking mass parameter Mq^ is equal for all members of an SU {2)l doublet, 
while the soft SUSY-breaking mass parameter M^^ can be different for scalar top and scalar 
bottom type quarks. Qq and denote the charge and isospin of q. Ag is the trilinear 
soft-breaking parameter, /i the Higgs superfield mixing parameter, tan/3 = V2/V1 denotes 
the ratio of the two vacuum expectation values in the Higgs sector (see Sec. 12.21) . Mz and 
Mw are the Z and W boson mass, respectively, and = cos 9^ = My//Mz with 6^ being 



the weak mixing angle, and s„ = y 1 — c^- The mass matrix can be diagonalized with the 
help of a unitary transformation Ug, 

Ti, = VrMM=M' ^ , U.= r^^" , fill 



-^-^^-^ y ml) ' \Uq.,, UqJ ' 

where the scalar quark masses, nig-^, nig^, will always be mass orderecjf] i.e. nig^ < rrig^, and 
are given by 



= i « + Ml) + ml + i/,^M|c2, 



T y[Ml -Ml + MlcpiPg - 2Qgsl)]' + AmlW ■ (12) 

For the parameter and the field renormalization of the quark/squark sector we follow the 
procedure described in Ref. [38]. The quark mass and the squark mass matrix are replaced 
by the renormalized mass and mass matrix, respectively, and their counterterms, 

rUg — )■ rUg + Snig , (13) 
M,- ^ M,- + 6Mg . (14) 

The mass matrix counterterm 6M.g is obtained by applying the renormalization procedure - 
replacement of the parameters by the renormalized ones and the corresponding counterterms 
- for each parameter and expanding with respect to the introduced counterterms, 

6Mg,, = 6 Ml + 2mg6mg - Mlc2p Qg 6sl + (/J - Q,^) (C2/3 5M| + M| ^2/3) , (15) 

(5Mg^2 = (^g ~ f^^) ^''TT'q + ^q{^^*q — f^Sfi — K 6fi) , (16) 

6Mg,, = 6Ml^ , (17) 
6Mg,, = 6 Ml + 2mg6mq + M|c2/3 Qq 5 si + Qgsl{c2p SM^ + M| 5c2p) (18) 

with K given in Eq. ffTOl) . 

Instead of starting out with the squark mass matrix in Eq. iQ the mass matrix in terms 
of the squark masses as given in Eq. (|TT1) can be used: 

U,M,ut^U,M,ut + U,5M,ut=(< II) ^ {^6Y* ^^9) 

where and are the counterterms of the squark masses squared. In the mass matrix 
in Eq. fllip as well as in the first term of the right hand side of Eq. (1191) the squark mixing 



Because of the mass ordering, &2 ~ is possible, which should be remembered when choosing a set of 
independent parameters - in our numerical examples, however, &2 is rather Br like. 
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parameter Yg vanishes as it should at tree-level because the unitary matrix Ug is chosen in 
that way to diagonalize the mass matrix M^. However, already at one-loop level, the squark 
mixing parameter Yg receives a nonvanishing counterterm 6Yg (which can be related to the 
counterterms of a mixing angle and a phase; see Ref. [39]). Using Eq. f|T9l) one can express 
SNlq by the counterterms Sm^^, and SYg. Especially for ^M^^^ and 5Mg^2 this yields 



SMg,, = \Uiy5mi^ + \UgJ'5mi - Ug,,U*gJYg - Ug,,U*gJY* (20) 
5Mg,, = Ul^Ug,,i6ml - 6ml) + Ul^Ug,,5Yg + Ug^.UiJYg* . (21) 



Equation fl2T|) can be used together with Eq. f|T6l) to express the counterterm for the trilinear 
top coupling 6 At and the counterterm for the bottom squark mixing parameter 6Yh by the 
other counterterms (see Eqs. fIFIl) and flS21) ). 

For the field renormahzation of the quark and the squark fields the following procedure 
is applied, 

uj.q^{l + l5Z!^)i0.q, (22) 
co+q-^il + ^SZ!^)co+q, (23) 
q^-,[l+lSZg].^q, . (24) 

[(5Zg]^_^. with i,j = 1,2 are the squark field renormahzation constants and 6Z^ and 6Z^ the 
field renormahzation constants for the left- and right-handed quark fields, respectively. 
Following this renormahzation procedure yields for the renormalized squark self-energies 

±g,, (/) = S,,, ip') + Up'- ml) [5Zg + 5Z|] - Sml , (25) 

tg,,{p') = S,,,(p2) + |(p2 - ml) [6Zg] + |(p2 - ml) [5Zt]^^ - 5Yg , (26) 

S,,, (/) = S,,, ip') + l{p'- ml) [6Zl] ^2 + i(/ - ml) [6Zg] - 6Yg* , (27) 

E,,, ip') = Eg,, {p') + l ip' - ml ) [6Zg + 5Z|] - 6ml . (28) 

The renormalized quark self-energy. Eg, can be decomposed into left /right-handed and scalar 
left /right-handed parts, Tig^^ and T,q^^^^, respectively, 

Eg(p) =^u;_Sj(p2) +^a;+Sj(p2) +a;_Sf (/) + (/) , (29) 

where the components are given by 

= Sj/^(p^) + l{6Zl/^ + 6Z^n , (30) 
Sf (P^) = ^'Ap') - "^i^Zl + 5Zf ) - 6mg , (31) 
(/) = Sf (p^) - ^(5Zf + 6Zi*) - 6mg . (32) 

It should be noted that ReEj^(p2) = {Ret^^{p^))* holds due to CVT mvariance. 

We now review our choice of renormahzation conditions where we follow the renormahza- 



tion scheme of Ref. [38] with our favored "mf,. Ah DR" scheme for the bottom quark/squark 
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part. However, we expand this scheme to include also external bottom quarks (which was 
not investigated in Ref. [M])- In this case we deviate from the "m^, Ajj DR" scheme and 
renormalize the bottom quark mass on-shell; see below. The problems found in Ref. [3B] 
with this scheme do not arise in the processes with external bottom quarks considered in 
this paper as no external bottom squarks occur, and the trilinear coupling Ab is only needed 
at leading order in these processes. 

The original parameters that we count as parameters of the top and bottom quark/squark 
sector are the soft SUSY-breaking mass parameters M^^, M^^, and M^^, the complex trilinear 
couplings At = \At\e^'^'^t, and Ai, = \Ab\e^'^'^t, and the Yukawa couplings yt and yt that can 
be chosen to be real (the Cabbibo-Kobayashi-Maskawa (CKM) matrix is set to unity in 
our calculation and generation mixing effects are neglected). Consequently, there are nine 
parameters to be defined in the top and bottom quark/squark sector. Instead of using the 
original parameters we choose the top squark masses m^-^ , m^-^ and one bottom squark mass, 

) ^^11 the quark masses rrit and mj, as independent input parameters]^ Also, in the 
scalar top quark sector a renormalization condition is chosen that fixes the counterterm 6Yt 
instead of 6 At. For the parameters of the top quark/squark sector we impose on-shell (OS) 
conditions while in the bottom quark/squark sector a mixed DR/OS scheme is employed: 

(i-iii) The two top squark masses and the one bottom squark mass are determined via on-shell 
conditions, 

ReE,-Jm|) = (2 = 1,2), (33) 
ReE. (mn = 0, (34) 



yielding 



Sml = ReEiJml) (^ = 1,2), (35) 
4 = f^eS^^^(m|; 



5m? =ReS. (m?) . (36) 



Re denotes the real part with respect to contributions from the loop integrals, but 
leaves the complex couplings unaffected]^ 

(iv) The top-quark mass is also defined on-shell, 

Re±t{p)t{p)l,^^,=0 , (37) 
yielding the one- loop counterterm 6mt, 

6mt = iRe {m, [S,^(m?) + Sf (m?)] + [Sf (m,^) + Sf^(m?)] } , (38) 
referring to the Lorentz decomposition of the self-energy St(p), see Eq. (l29l) . 



^ It should be noted that in the case &2 ~ fni^ cannot be chosen as an independent parameter and 
the renormahzation scheme has to be switched to one with m^^^ as an input parameter. 

^ It should be noted that we impose later an extra renormalization condition concerning the &i-mass to 
solve infrared problems, see Eq. (^5]) below. 
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(v) For the bottom quark mass we use two different definitions depending on whether in 
the considered decay channel a bottom quark appears as an external particle or not: 

(a) In the case that no external bottom quarks are involved the bottom quark mass 
is defined as DR-mass with the corresponding counterterm, 

^mp = |Re {m, [S^(m^) + Sf (m^)]^^^ + [Sf K) + Sf^(m,^)]^^J . (39) 

(b) If bottom quarks appear as external particles then we define the bottom quark 
mass on-shell, 

ReSfe(p)6(p)|^,^^, = 0, (40) 

to ensure the on-shell properties of the external particles which yields the following 
counterterm: 

Smf^ = |Re {m, [S^(m^) + (m,^)] + [Sf (m^) + Sf«(m^)] } . (41) 

To have consistent input in all the decay channels we calculate the on-shell bottom 
quark mass m^^ starting from the DR-mass, 

= mf^ + 6mf^ - 6m^^ . (42) 

The value of is obtained as described in Eq. fl65|) . 

It should be noted that the problems found in Ref. [3H] with an on-shell renormalization 
condition for (leading, e.g., to unphysically large contributions to 6Ai,) do not 
occur as long as no external scalar bottom quarks appear at the same time and the 
parameter Ai, is only needed at leading order. For example, the proposed scheme would 
presumably fail in the process bb — j- biij, which, however, is beyond the scope of our 
paper. 

(vi,vii) The complex counterterm of the non-diagonal entry of Eq. (ITU]) , which corresponds to 
two separate conditions, is fixed as 



6Y, = iRe {^a^D + S,-,,(m|)} . (43) 

(viii,ix) In the scalar bottom quark sector the trilinear coupling is defined as a DR parameter 
with the counterterm, 



SAb = — 
nib L 



+ -2 Ul^M.., (ReS^^^(m[)Uiv + ReS^^^(m|)Uiv 
-\{Ab- fi* tan/3) Re{m, [S^(m,^) + (m^)]^.^ 
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+ SfJ'*\div tan f3 + fi* (5tan/3 , (44) 

which also counts for two separate renormahzation conditions as is a complex pa- 
rameter. The divergent parts of 5/i and (5tan/3 can be extracted from the Eq. (11811) 
and Eq. (I120cp . respectively. 

With these renormahzation conditions all independent parameters in the top and bottom 
quark/squark sector are defined. The dependent parameters can be expressed in terms 
of those independent ones and the same applies for the corresponding counterterms. With 
respect to this renormahzation scheme the (one- loop corrected) on-shell bi mass, m?^, differs 

from the mass parameter . As an external particle bi should fulfill the on-shell properties, 
which in turn requires that it should have the mass 



Kr = K) + (KD - ReS,^^(m[) (45) 



where Sm^^^' is the dependent mass counterterno that results from imposing only the renor- 
mahzation conditions (i) - (ix). On the other hand, using m^^ for the internal bi squarks and 
^os Jqj, ^Yie external bi squarks (which would formally be correct with respect to the consid- 
ered loop-order) leads to nonvanishing infrared (IR) singularities (for details see Ref. [38]). 
To circumvent this problem we impose a further OS renormahzation condition 

Sml = ReE-,Jml) . (46) 

As now all the squark masses within one generation are renormalized as on-shell an explicit 
restoration of the SU{2) relation is needed. This is performed in requiring that the left- 
handed (bare) soft SUSY-breaking mass parameter is the same in the bottom as in the top 
squark sector at the one-loop level. 

Ml (b) + 5 Ml (b) = Ml (i) + 6 Ml (i) . (47) 

More precisely, we define (see also Refs. [TSllT^I^ ) 

Ml (6) = Ml (i) + 5 Ml (i) - 5 Ml (b) (48) 

with 

6Ml{q) = \U,J'5ml + \U^J'5ml - U.^.U^JY, - U,,,U;jY* - 2m,6m, 

+ Ml Q, 6 si - ( Jj -Q,sl) {C2p 5 Ml + M| dc^p) . (49) 

where (5M|^(g) is derived with the help of Eqs. f lT5|l and f l20|) . Now M|^(6) is used in the 
scalar bottom mass matrix instead of the parameter Ml in Eq. (jH]) when calculating the 
values of m^^ and m-^^. However, with this procedure, also the mass of the 62 squark is 



It can be found, up to 0{as) in Eq. (8.43) of Ref. 
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shifted, which contradicts our choice of independent parameters. To keep this choice, also 
the right-handed soft SUSY-breaking mass parameter M^^ receives a shiftEl: 

Ml = . mg|A*-/itan/3p ml - c^s Qb si + m? (50) 

Ml{b) + ml + Mlc,,{n-Q,sl)-ml " ^'^"^'^ 

Taking into account this shift in M^^, up to one- loop ordej^, the resulting mass parameter 
mi is the same as the on-shell mass Eq. fHSl) . 

In the top and bottom quark/squark sector the counterterm for the trilinear top coupling 
SAf and the counterterm SYt, are given as a combination of the independent parameters that 
can be derived from the relation of Eqs. ( IT6|) and (12T]) . 

SA = ^[Ui^^Ul^{6ml-6ml) + Ui^^Ul^6Y: + Ul^^^^^^ 

I ill 

+ (5/i* cot /3 - /i* cot^ /3 5tan/3 ) (51) 



and 

5Y. 



Ur m (6ml -6ml 

Oil 621 V bi 62 



+ mh {Ui^Ul^^ {6AI - ji 5tan/3 - tan /3 5/i) - Ui^^U^^^ {6Ab - fi*6taii(3 - tan f3 6fi*) ^ 
+ (t/^^^f/4^(A*-^tan/3)-f/^^^[/^^^(A-^*tan/3)) 6m^] , (52) 

where 5tan/3 and 6fi will be defined within the Higgs/gauge sector in Sec. 12.2^ Eq. f ll20cp 
and the chargino/neutralino sector in Sec. 12. 3[ Eq. fllSip . respectively. 

Now, the parameter renormalization for the top and bottom quark/squark sector is ac- 
complished but the field renormalization still has to be done. We determine the Z factors of 
the quark and squark fields in the OS scheme. In the quark sector we have 

ReS,(p)g(p)|^2=„. = 0, (53) 

lim ^ '^ \ ip) = 0, 54 

where these two equations determine not only the quark mass counterterms (see Eqs. fl371) 
and fHOj) ) and the real part of the Z factors but also the difference of the imaginary parts 
of the quark Z factors, Im 6Zg — \m.6Z^ (analogously to the chargino/neutralino case in 
Sec. 12. 3p . This leaves us the freedom to impose additionally 

\x^6Z^ = -\^6Z^ . (55) 



^ If the mass of the h\ squark is chosen as independent mass as hi « 6l then the shift of M^^ has to be 
performed with respect to mg^ . 

^In the case of a pure OS scheme (see e.g. [42ll43] for the rMSSM) the shifts Eqs. gH]) and dSO]) rcsuh in 
a mass parameter m~^^^ which is exactly the same as in Eq. (|45p . This constitutes an important consistency 
check of these two diiferent methods. 
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With these equations we find 



Re6Z^/^ = -Re 



+ m„ 



R' 



K) + sf V^) 



Im^Z^/^ = ±^Re jSf^K) - Sf K)| = ±^ Im (ReS 



For the scalar quarks we demand 



(z = l,2) 



ReSq,2(m?J = ReSg2,(m? 



92iV"''gi; 



yielding 



Re[5Z,]^^ = -ReS^Jp2 



1,2) 



12 



21 



+2 



-2 



ReS,-,,(m|)-^y, 
(m? - m? ) 



ReS,-,,(m?J-^y; 
(m? - m? ) 



with S'(p^) = ^§^^5 (1 = {^7^}; and q = {t,b}. It should be noted that the on-shell 
conditions leave the imaginary part of [5Zq] .. undefined; it can be (implicitly) set to zero as 
it does not contain any divergences, 



(56) 
(57) 

(58) 
(59) 

(60) 



(61) 



Im[5Z,]^ 







1,2) 



(62) 



The input parameters in the b/b sector have to correspond to the chosen renormalization. 
We start by defining the bottom mass, where the experimental input is the SM MS mass 



mf^{mb) = 4.2 GeV 



(63) 



The value of m^^{fifj) (at the renormalization scale = m^J is calculated from mf^^{mh) 
at the three- loop level following the prescription given in Ref. |15]. 
The "on-shell" mass is connected to the MS mass via 



mi 



1 + ffll l + 21n^^ 

3 ,„MS(^^) 



(64) 



The DR bottom quark mass at the scale /i/j is calculated iteratively from | ^ H3 | H6] 



m 



DR 



^os 


l + Afc 


1 + Sm^^ - Smf^ 




1 + A, 





(65) 
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with an accuracy of 1 1 — (mf^) ^"•'^ / (mf^) | < 10"^ reached in the nth step of the iteration 
where 6m^^ and 6m^^ are given in Eqs. f l39l) and fHT]) . respectively. 

The quantity [^61418] resums the C((a<itan/3)") and tan/3)") terms and is given 

by 

A. = tan/3M3>*/KVm|,m^) + tan/3 /(m|, m|, |/ip) (66) 

with 

, , ab\n{b/a) + ac\ia.{a/c) + bc\n{c/b) , , 

i[a,b,c) = - r- — — - . (67j 

[a — c)[c — b)[b — a) 

Here at is defined in terms of the top Yukawa coupling yt{rnt) = \/2mt{mt)/v as at{mt) = 



y'iimt) / {A-k) with v = 1/vV2Gf = 246.218 GeV and mt{mt) ^ mt/{l - j^^atimt) + 
-^as{mt)). Setting in the evaluation of Af, the scale to nit was shown to yield in general a 
more stable result [19] as long as two-loop corrections to A;, are not includedl^. M3 is the soft 
SUSY-breaking parameter for the gluinos; see below. We have neglected any CKM mixing 
of the quarks. 

2.1.2 The gluino sector 

The gluinos appear as external particles, for instance, in the decay ^2 tg. Therefore, a 
renormalization procedure for the gluino field and the corresponding parameters is necessary. 
The Fourier transformed Lagrangian bilinear in the gluino fields is given by 

9org = 9ovgf ^org + 9ovgf W+ ^org " ^org ^3 ^- 9org " ^org ^3 ^+ ^org (68) 

with M3 being the soft-breaking gluino mass parameter, which is in general complex, 

Ms = \M^\e''^'^ . (69) 
The gluino field g'^^^ can be redefined using a phase transformation 

uJ±r = e^'^c^i^^rg (70) 

such that the gluino phase ^Pg appears only in the gluino couplings, but not in the mass term 
with the gluino mass nig = {M^l. 

The renormalization is performed as follows [50] : 



M3 -> Ms + 5Ms = Ms + (5|Ms|e*'^s + iM^Sip-g , 
u^r ^ (1 + ISZ-g) oo^r , 

r ^ (1 + uj^r ■ (71) 



^ It should be noted that in Ref. [49| a different scale has been advocated due to the emphasis on the 
two-loop contributions presented in this paper. The plots, however, show that mt is a good scale choice if 
only one-loop corrections are included. 
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Note that, analogous to the mixing matrix, only the gluino phase appearing in M3 is renor- 
malized but not the one appearing due to the redefinition of the gluino field. 
The renormalized gluino self-energies read 

= ""(P') + WZ, + 5Zl) , (72) 
Ef (/) = Sf _ 5Z-^ - 6M,e-'^-^ , (73) 
Ef (p2) = Ef (p2) - SZl - 5M;e'^^ . (74) 

We choose OS renormalization conditions for the gluino, 

Ret-,{p)r{p)\ = (75) 

p2^-m.? p — m~ 

with S(p) defined according to Eq. (12^ . Because of the Majorana nature of the gluino this 
leads to three independent conditions, yielding - with E'(m^) = ^§^-^1^2-^2 ^ 

5IM3I = ^R^ {m, [Ef (m|) + Ef (m|)] + [Ef (m^) + Ef (m|)] } , (77) 
Re5Z~, = -Re{Ef (m|) + m| [Ef (m^) + Ef (m|)] + m~, [Ef' (m|) + Ef' (m|)] } , (78) 
lm5Z, = ^Re {Ef (m|) - Ef (m^)} - 6^, . (79) 

We have then chosen 6(pg = 0, which is similar to the quark case: There, we use a real 
Yukawa coupling due to the possibility of redefining the quark fields and have a complex 
Z factor at one-loop order that keeps the Yukawa coupling real also at one-loop order; in 
contrast the gluino phase still appears in the Lagrangian after the redefinition of the fields 
but this phase factor can be considered as a "transformation matrix" and does not obtain 
a counterterm. Note that in the chargino/neutralino sector we keep the diagonal Z factors 
real and have complex parameters. 

2.1.3 The strong coupling constant 

The strong coupling constant is renormalized as 

as -> Za.as = (1 + SZaJas ■ (80) 
For Za^ we are using a DR renormalization condition yielding 

5Z^^ = -lEl'^{p%,^^^^^^, (81) 

where E^g denotes the derivative of the transverse part of the gluon self-energy. 

The decoupling of the heavy particles and the running is taken into account in the 
definition of as', starting point is [H] 

ap'(^)(Mz) = 0.1176 , (82) 
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where the running of a^^'^"'^\fifi) can be found in Ref. [H]. From the MS value the DR 
value is obtained at the two-loop level via 



a 



DR,(n/) 



a 



MS,(n/) 



MS,(n.^) 
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For fifi > rrit we have Uf = Q. Within the MSSM, at one-loop level, as reads 



a 



MSSM 



a 



DR,(6) 



R 



a. 



DR,(6) 



TT 



In^ + lni^ 



^3) 



^4) 



with Mg being defined as the geometric average of all squark masses, Mg = Ilg{mqj^mq^)'i^ . 
The log terms originates from the decoupling of the SQCD particles from the running of 
at lower scales fiR < /idee. = ^i2- -^^^ simplification we have chosen m^^ (representing the 
energy scale of the considered decays and as a typical SUSY scale) also as decoupling scale. 



2.2 The Higgs and gauge boson sector of the cMSSM 

The MSSM Higgs potential Vh, 

+ \{gl + gDmn,, - H;,H^,f + \gl\H*M' (85) 

with = {1,2} and e^^ = 1, contains both the U (1) and SU (2)^ gauge coupling constants 
Qi and g2, respectively, which are considered to be part of the gauge boson sector as well as 
the soft SUSY-breaking parameters mi2, rh\, and (with m\ = mf + = rfi^ + |/ip), 

which are part of the Higgs sector. For this reason we do not separate those two sectors but 
treat them within one section. The Hij with {i,j} = {1,2} are the components of the two 
Higgs doublets that can be decomposed in the following way: 



Hl2 J \ — 

Besides the vacuum expectation values V\ and U2, in Eq. (186|1 a possible new phase ^ between 
the two Higgs doublets is introduced. 

In total, the Higgs and gauge boson sector contains 7 real parameters: g^-, ml, vi, 
V2, and ^ and one complex one mi2. /i is defined within the chargino/neutralino sector; see 
Sec. 12.31 With the help of a Peccei-Quinn transformation |52] /i and can be redefined [53] 
such that the complex phase of vanishes. 

The part of the Fourier transformed Lagrangian that is linear or bilinear in the massive 
gauge boson and Higgs boson fields is the following: 

^ffiggs = 'Pi + 02 + Xi + T^2 X2 



Corresponding to the convention used in FeynArts/FormCalc, wc exchanged in the charged part the 
positive Higgs fields with the negative ones, which is in contrast to [37]. As we keep the definition of the 
matrix M^i^i used in |37| the transposed matrix wih appear in the expression for M^ljfgj-; see below. 
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(01,02,Xl,X2) (p'll-M, 



<l'4>xx) 



92 



fii/ 



02 

Xi 

\X2j 



w;:[{p' - M^)g^^ - (1 - 

^2 



p^Mw{W+G- + G+W-) - ^wM^G^G- , 



^7) 



where we have used the decomposition of Eq. (186|) . The coefficients of the hnear terms are 



called tadpoles with Ts 



{1,2} being the tadpole parameters. M,^,^^^ and M,^ 



are the Higgs mass matrices. The terms containing the neutral and the charged Goldstone 
boson fields, G and G^, and the gauge parameters, and respectively, are coming from 
the gauge-fixing part of the Lagrangian. The interaction fields can be transformed into mass 
eigenstates (see Refs. [37lll0],[5l] for details, especially on notation) with the help of the real 
and orthogonal transformation matrices Un(o) and Uc(o) , 



(h\ 






H 

A 


= Un(0) 


4>2 

Xi 






\X2j 




) = Uc(0) 


(1") 



With Mf^\^ 



and {Th,TH,TA,TG) — (T^^^, T^j, T^^, T^jUj^^g^ 



with M 



diag 
H±G± 



u, 



c(0)iVl^±^i 



(89) 
(90) 



where the diagonal elements of M^'^^ig ^'h±g± tree- level masses denoted as rrih, 

TTiH, ^A, and Mh±, mG±, respectively. It should be noted that the tadpole parameter 
Tg can be expressed by tan/3, the entries of Un(o) and Ta- 

Throughout our calculation we use the 't Hooft-Feynman gauge, C,z = = 1- Concern- 
ing the renormalization procedure, we follow the usual approach where the gauge-fixing terms 
do not receive a net contribution from the renormalization transformations. Accordingly, no 
counterterms as given below arise from the gauge-fixing terms. 



We replace the 8 original parameters vi, V2, gi, g2, ^ 



2 

12 5 



;^2 ^2 
2i 



and ^ by the Z and 



W boson mass Mz, Mw, the electric charge e, tan (3, the mass of the charged Higgs boson 
Mh±i and the tadpole parameters T/i, Th, and Ta (where we have chosen to be real, 
which is always possible with the help of a Peccei-Quinn transformation). Details about this 
replacement can be found in Ref. [37] . 

The minimization of the Higgs potential in lowest order leads to the requirement that 
the tadpole coefficients Tf^h H,A} in Eq- (IHTl) must vanish (the tadpole coefficient Tg vanishes 
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automatically if = holds, as Tq can be written in terms of T^). In particular, the 
condition T^i = implies that the complex phase ^ has to vanish, see e.g. Refs. [HZj, so that 
the Higgs sector in lowest order is CV conserving. 

In order to derive the counterterms entering the one-loop corrections to the Higgs boson 
masses and effective couplings, we renormalize the parameters appearing in the linear and 
bilinear terms of the Higgs potential^. 



e 



Mi 



;i + 5^e)e , 
tan /3 — )■ tan /3 + 5tan/3 



(91) 



Th^Th + 6Th , 
Ta^Ta + 6Ta . 

It is important that according to our renormalization procedure the renormalization prescrip- 
tion has to be applied before the transformation into the mass eigenstates, also for tan /3, i.e. 
a P appearing from the transformation to the fields A, G, H^, and does not obtain a 
counterterm. For the counterterms arising from the mass matrices we use the definitions 



(92) 
(93) 



with 



6M 



hHAG 



Un(0) 5M, 



/ 5ml 




^^Ig\ 


5mljj 








^^HA StTT'A 








5ml I 


( 5m\± 


^^G-H+\ 






+ 6ml± J 





(94) 



(95) 



where (5M, 



'P'f'XX 



and (5M^±^± denote the counterterm mass matrices that are obtained when 



replacing the parameters in M,^,^^^ and M^^i^i by the renormalized ones and their coun- 
terterms using Eq. fpTl) . expanding and taking the first-order expressions and applying the 
zeroth order relation T!^h,H,A} = 0; see Ref. [37] . 

As mentioned above, in contrast to what is often done in the MSSM with real parameters, 
we use Mh± as an independent input parameter. The counterterm 6m\ in the formulas above 
is therefore a dependent quantity, which has to be expressed in terms of using 



6m\ 



Sm^ 



(96) 



For the field renormalization we choose to give each Higgs doublet one single renormal- 
ization constant, 

n,^il + \5Zh,)V-i , ^2 ^ (1 + \SZn,)U2 . (97) 

-'^^ It should be noted that in Ref. [37] a shghtly different renormahzation prescription for tan/3 had been 
introduced, tan /3 — !• tan /3(1 + (Stan/3 such that (5tan/3 = tan /3 (5tan/3 
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In the mass eigenstate basis, the field renormahzation matrices read 



/h\ 




/h\ 


H 




H 








A 




A 









1 

+ 2 



/ SZhh SZhH SZhA SZhG\ 

SZhH SZhh SZha SZhg 

5ZhA 5Zha 5Zaa SZag 

\5ZhG 5Zhg 5Zag SZggJ 



H 

A 

\GJ 



(98a) 



and 




(98b) 



5Zg-g+ cire not independent but can be derived 



where, according to Eq. ( 197|) . 5Zhh, ■ 

via Un(o) diag(5Z-Hi, 5^-H2, ^^Wi, ^^Wa) ^\(o) and Un(o) diag(5Z^j, ^ yielding the 



following expressions for the field renormalization constants in Eq. (!98lll^ : 



8Zhh 


= sin^a + cos^a SZ-}^^ , 


(99a) 


5Zaa 


= sin^/3 SZ-}i-^ + cos^/3 SZ-}^^ , 


(99b) 


5ZhH 


= sin a cos a {SZ-}^^ — SZy^J , 


(99c) 


SZag 


= sin/3 cos f3{6Z-H2 - SZ-^J , 


(99d) 


5Zhh 


= cos^a + sin^a SZ-}^^ , 


(99e) 


8Zgg 


= cos^/3 + sin^/3 SZ-}i^ , 


(99f) 


SZh~h+ 


= SZaa , 


(99g) 


SZh-g+ 


= SZg-h+ = SZag , 


(99h) 


SZg-g+ 


= SZgg ■ 


(99i) 



For the field renormalization constants of the CP-violating self-energies it follows, 

SZhA = SZhG = SZha = SZhg = , 



(100) 



which is related to the fact that the Higgs potential is CV conserving in lowest order. 

In the case of a decay to a neutral or charged Higgs boson some self-energy transitions 
on the external Higgs leg have to be taken into account. In order to define the various coun- 
terterm contributions, we list here the respective renormalized self-energies (taking already 
Eq. fllOOl) into account). For a scalar-vector self-energy we use Sgy(p'^) = p'^Ssv(p^), where 
p'^ is the momentum of the incoming scalar or vector particle. Then we have the following 
renormalized self-energies: 



-•HG 



(101) 
(102) 



It should be noted that Eq. ( |99g[ ) is sufficient to yield an UV finite result for decays involving a charged 
Higgs boson. To obtain also an IR finite result and to ensure the on-shell properties of the outgoing charged 
Higgs boson, another Z factor has to be taken into account; see Eq. (|134l) below. 
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^ag{p ) 


v f 2\ , <: ry / 2 1 2 \ r2 

= SagIP ) + (>Zag[P - - om^G , 


^H-G+[P ) 


= j:h-g+{p > + oZh-g+[p - i^H±) - 


^hz{p^) 


= ^hzip") , 


f^Hzip^) 


= ^Hzip') , 


tAzip') 


= T.Azip'^) - Sm\z , 


^H-W+ip"^) 





(103) 
(104) 
(105) 
(106) 
(107) 
(108) 



with the mass counterterms expressed as [37] 



SmlG = " ""f " , (109) 

'^^ 2Mzs„c^ ^ ' ^ ^ 

Sml^ = -'-^^^^^^5TA, (110) 

^ 9 e (sin(/3 — a) — cos(/3 — a) 5T/i) 9^/,^9 , ,9 x r ^ /...x 
5m^g = ^- '—^ - cos^(3 (M^± - M^) <5tan/3 , (111) 

2 e (sin(/3-a)(5TH-cos(/3-a)(5Tft + z5TA) 20.^2 o n1o^ 

6mH-G+ = cos /3 M^± (5tan/3 , (112) 

(5m^2 = +iM^ cos /3 [cos /3 5tan/3 + | sin /3((5Z-^2 ~<^^>^i)] j (113) 
5m^-^+ = — Mvi/cos/3 [cos/3 5tan/3 + | sin /3((5Z-^2 — 5^>^J] . (114) 

a denotes the angle that diagonahzes the CP-even Higgs boson mass matrix at the tree-level. 
It should be noted that according to Eq. (157|) , no mixing of the CP-even Higgs fields and the 
Z boson fields occurs at tree-level. Consequently, there are no counterterm contributions to 
this mixing at one-loop level; see Eqs. (11051) and (11061)1^ 

In the following we list our renormalization conditions and the resulting counterterms: 

(i-iii) We impose on-shell renormalization conditions for the masses of the SM gauge bosons 
and the charged Higgs boson, 

R^EL(M|) = 0, R^S^i^(M^) = , ReS^^+^,-(M^±) = . (115) 

The gauge boson self-energies S-^ are the transverse parts of the full self-energies. 
Eq. flll5p yields for the mass counterterms, 

5M| = ReSi^(M|) , <5M^ = ReE^^(M^) , Sm^^ = ReEH+H-{Ml±) . (116) 



(iv) The electric charge is defined via the standard on-shell conditions requiring that no 
corrections occur to the electron-positron-photon vertex with on-shell external particles 

The other renormalized Higgs boson self-cncrgics, T,hh, ^hh, ^hh, ^aAi ^gg, ^hA, ^ha, and T,h-h+ 
and the corresponding mass matrix counterterm contributions, Sm^, <5m|^, Smf^j^, Stuq, 5m\j^, ^'^'hai 
are not expHcitly needed in our calculation (employing FeynHiggs, these self-energy contributions are au- 
tomatically taken care of). They are given by Eqs. (64) and (53) of Ref. [37], respectively, considering 
5tan/3 133 = Stai^P / tan p. 
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at zero photon momentum. This yields bZ^ expressed through the photon and photon- 
Z self-energies, 

(v-vii) The tadpole parameters are renormalized such that the complete one-loop tadpole 
contributions vanish, 

^Sk^} + '^^{^.^,^} = , (118) 

where t'^^^^j^^ denote the contributions coming from the genuine one-loop tadpole 
graphs. These conditions lead to 

= -Tt'^ , STh = -T^H^ , 5Ta = -T^A^ . (119) 

(viii) The last parameter that has to be defined is tan /3. We do that together with the Higgs 
boson field renormalization constants 5Z-^^ and SZy^^. 

A convenient choice that avoids large (and unphysically) higher-order corrections in 
the (c)MSSM Higgs sector is a DR renormalization of SZy^-^^, SZy^^, and (5tan/3 [55] . 



SZn, ^ SZ^f = -ReSU(0)L^o ,,, , (120a) 

6Zn, ^ 5Z^=-ReKhiO)l=o4i. ' (120b) 
5tan/3 = 5tan/3 °^ = 1 tan/3 (5^^^ - ^^wJ , (120c) 
i.e. the counterterms in Eqs. f ll20p contribute only via divergent part^, and the finite 



result depends on the renormalization scale fiR. For the setting of fiR see Sec. 14.11 

The DR renormalization of the parameter tan /3, which is manifestly process indepen- 
dent, is convenient since there is no obvious relation of this parameter to a specific 
physical observable that would favor a particular on-shell definition. Furthermore, the 
DR renormalization of tan/3 has been shown to yield stable numerical results [551457] . 
This scheme is also gauge independent at the one-loop level within the class of 
gauges 

Finally, the field renormalization constants of the gauge bosons have to be determined. 
Applying an on-shell condition for the gauge boson fields, the field renormalization constants 
can be derived as 

SZzz = -ReSFz(Ml) , 6Zww = -ReJ^^wiM^) ■ (121) 
In other sectors, one may need the following counterterms expressed by independent ones: 

5 sin /3 = cos'^/3 (5tan/3 , 5 cos /3 = — sin /3 cos^/3 (5tan/3 . (123) 

■^^ The divergences in Eqs. (|120ap . (jl20b|) arc momentum independent. 
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We have checked (at the one-loop level) that the following Slavnov- Taylor identities [581 - 
hold: 

^kcip') - '-^^hzip') = , (124) 



Shg(p') - -^^Hzip') = , (125) 



^2 



Sag(p') - '-^^Azip') + ip' - ml) Up') = , (126) 



2ip^<^ , 2^ 

,2\ P ^ /2\ , / 2 71^2 \ _f / 2 



Sgg(p^) - -^^Gzip") - ^S|^(p2) = , (127) 

2 



^h-gAp') - ^SH-iy+(p') + ip' - M^±)/±(p^) = , (128) 



^g-gAp') - ^^g-wAp') - -^^wwip') = , (129) 



where S"^ denotes the longitudinal part of the self-energy and 

/o(p') = "2 2 sin(/3 - «) cos(/3 - a) [B,{p\ ml M|) - Bo{p', m\, M^] , (130) 

f±{p') = -TT^ sin(/3 - a) cos(/3 - a) [5o(p^ M^) - 5o(p^ m^, M^)] . (131) 



The definition for the Bq function can be found in Ref. [61]. The Slavnov- Taylor identities 
also hold for the unrenormalized self-energies (where the tadpole contributions must not be 
neglected) . 

The Higgs boson field renormalization constants are necessary to render the one-loop 
calculations of partial decay widths with external Higgs bosons UV finite. The DR scheme 
for the field renormalization constants is used in the calculation of the Higgs masses within 
FeynHiggs in order to avoid the possible occurrence of unphysical threshold effects. As the 
results of FeynHiggs are used within the numerical evaluation in Sec. HI it is appropriate and 
consistent to follow the same renormalization procedure. As always, Higgs bosons appearing 
as external particles in a physical process have to obey proper on-shell conditions. A vertex 
with an external on-shell Higgs boson hn {n = 1,2,3), Th^, is obtained from the tree-level 
vertices Th, Th, and Ta via the complex matrix Z [37] . 

= [ZUTh + [ZUTh + [ZUTa + ... , (132) 

where the ellipsis represents contributions from the mixing with the Goldstone boson and 
the Z boson; see Sec. [3l It should be noted that the transformation with Z is not a unitary 
transformation; see Ref. [37] for details. 

Also the charged Higgs boson appearing as an external particle in a t decay has to obey 
the proper on-shell conditions. The corrections to the charged Higgs boson propagator lead 
to an extra Z factor, 



ReS'^-^+(/)| " . (133) 

' H± 
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Expanding both sides of Eq. fll33p up to one- loop order and using Zh-h+ = 1 + SZh-h+, 
leads to 

6Zh~h+ = -ReS'^_^+(/)| , = -ReE'^_^+(M^±) - 6Zh-h+ • (134) 
Analogous to the procedure for the neutral Higgs bosons; see Ref. [37] . 

^Jzh-h+ = 1 + IsZh-h^ (135) 

has to be applied to a process with an external charged Higgs boson. Within the presented 
calculations, for the charged Higgs bosons, we include contributions from sj Zh-h+ strictly 

at the one-loop level i.e. the correction coming from ^5Zh-h+ of Eq- (I134p multiplied by 
the corresponding tree-level vertex contribution]^ As for the neutral Higgs bosons, there 
are contributions from the mixing with the Goldstone boson and the W boson, which we 
deal with separately calculating the mixing explicitly and strictly at one-loop order. The 
Z factor Zfj-H+ is UV finite by definition. However, it contains IR divergences that cancel 
with (IR divergent) soft photon contributions from the one-loop diagrams; see Sec. [31 



2.3 The chargino/neutralino sector of the cMSSM 

The chargino/neutralino sector contains two soft SUSY-breaking gaugino mass parameters 
Ml and M2 corresponding to the bino and the wino fields, respectively, as well as the Higgs 
superfield mixing parameter /i, which, in general, can be complex^ The gauge boson masses 
and tan /3 that also appear in this sector have already been defined within the context of the 
Higgs and gauge boson sector; see Sec. 12.21 For our calculation we also need to renormalize 
the chargino and neutralino fields. 

The starting point for the renormalization procedure of the chargino/neutralino sector is 
the part of the Fourier transformed MSSM Lagrangian that is bilinear in the chargino and 
neutralino fields, 

4-, x« = + ^+*^^ - ^ [V^X^Ut],, co.xj - ^ [UX*VT],,. co+xj 

+ ^ co.xl + ^+Xl - ^ [N* YNt],, - 1 [NY*N^],., , (136) 

already expressed in terms of the chargino and neutralino mass eigenstates xi and xt^ 
respectively, and i,j = 1,2 and k,l = 1,2,3,4. The mass eigenstates can be determined 
via unitary transformations where the corresponding matrices diagonalize the chargino and 
neutralino mass matrix, X and Y, respectively. 

In the chargino case, two 2x2 matrices U and V are necessary for the diagonalization 

In our calculational set up we add Eq. (|134p to Eq. p9gP . 

Often, Af2 is chosen to be real, which is possible without loss of generality as not all the possible phases 
of the MSSM Lagrangian are physical and there is a certain freedom of choice; see the discussion in Sec. 14.31 
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of the chargino mass matrijo X 



V* u"f 







with X 



M2 

y/2 cos /3Mw 



V2 sin (3 Mw 



(137) 



where Mv- is the diagonal mass matrix with the chargino masses m -± , m-± as entries, 

^ ^ Xi X2 

which are determined as the (real and positive) singular values of X. The singular value 
decomposition of X also yields results for U and V. Using the transformation matrices 
U and V, the interaction Higgsino and Wino spinors iff, H2 and W^, which are two 
component Weyl spinors, can be transformed into the mass eigenstates 



with = U, 



w- 



and ipf = Vi 



w- 



(138) 



where the ith mass eigenstate can be expressed in terms of either the Weyl spinors ip^ and 
il^f" or the Dirac spinor . 

In the neutralino the neutralino mass matrix Y is symmetric, one 4x4 matrix 

is sufficient for the diagonalization 



M^o = N* Y = diag(m^o , m^o , m^o^ , m^^o j 



with 



Ml 

M2 

—Mz Sw cos 13 Mz Cw cos /3 

y M^SwSin/3 -M^c„ sin /3 



— Mz cos [3 Mzs^ sin [3 \ 

Mz Cw cos /3 — Mz Cw sin [3 

-II 

-/i 



(139) 



(140) 



The unitary 4x4 matrix N and the physical neutralino masses m^o (fc = 1, 2, 3, 4) result from 



a numerical Takagi factorization |62] of Y. Starting from the original bino/wino/Higgsino 
basis, the mass eigenstates can be determined with the help of the transformation matrix 
N, 



with 



kl 



:i4i) 



where denotes the two-component Weyl spinor and xt four-component Majorana 
spinor of the fcth neutralino field. 

Concerning the renormalization we follow the prescription of Ref. [50]. The following 
replacements of the parameters and the fields are performed according to the multiplicative 
renormalization procedure: 



Ml Ml + SMi 



(142) 



Corresponding to the convention used in FeynArts/FormCalc, we express the chargino part in terms of 
negative chargino fields, which is in contrast to |50j . As we keep the commonly used definition of the matrix 
X the transposed matrix appears in the expression for M^- . 
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Ma Ma + SM2 



J 



^+Xi 



kl 



(^,J = 1,2) , 
(^,J = 1,2) , 
(A;,/ = 1,2,3,4) , 
(A;,/ = 1,2,3,4) . 



(143) 
(144) 

(145) 
(146) 
(147) 
(148) 



It should be noted that the parameter counterterms are complex counterterms that each 
need two renormalization conditions to be fixed. The transformation matrices are not renor- 
malized, so that, using the notation of replacing a matrix by its renormalized matrix and a 
counterterm matrix 



X ^ X + (5X 
Y ^ Y + (5Y 



(149) 
(150) 



with 
5X 

5Y 



( 



5M2 
V2S{Mwcos(3) 





V2 5{Mwsm(]) 
Sfj, 



5M2 



-6{Mzs„ COS /3) 6{Mzc^ cos (3) 
\ (5(MzSwSin/3) - 5 {Mzc^ sin (5) 



-S{Mzs^ cos (3) 6{Mzs^ sin (3) \ 

S{Mzc^ cos (3) -6{Mzc„ sin (3) 


-Sji / 



the replacements of the matrices M-j^- and M-j^o can be expressed as 

M^o ^ M^o + 6M^o = M^o + N*5YN"f . 
Now the renormalized self-energies are given by 



(151) 



(152) 



(153) 
(154) 



E|«(/) 



Sfo(/) 



J kl 



J kl 



Sf(p^) 
Ef(p^) 



5Zj_ + 5Zjt 



i5Zf!M,-+lM,-5Z^-+5M, 



Sfo(/) 



kl 



kl 



2 X X 



6Z^o + 6Zl 



dZ^o + 5Z lo 



2 X X ' X 



(155) 
(156) 
(157) 
(158) 

(159) 
(160) 
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kl 



kl 



kl 



kl 



J kl 



kl 



(161) 
(162) 



Instead of choosing the three complex parameters Mi, M2, and /i as independent param- 
eters, we impose on-shell conditions for the two chargino masses and the mass of the hghtest 
neutrahno and extract the expressions for the counterterms of Mi, M2, and accordingly. 
In a recent analysis [63] it was emphasized that in the case of the renormalization of two 
chargino and one neutralino mass, always the most binolike neutralino has to be renormalized 
in order to find a numerically stable result. Also, in Ref. [61] the problem of large unphysical 
contributions due to a non-binolike lightest neutralino is discussed. In our numerical set up, 
see Sec. HI the lightest neutralino is always rather binolike. On the other hand, it would be 
trivial to change our prescription from the lightest neutralino to any other neutralino. In 
Ref. [63] it was also suggested that the numerically most stable result is obtained via the 
renormalization of one chargino and two neutralinos. However, in our approach, this choice 
leads to IR divergences, since the chargino mass changes (from the tree-level mass to the 
one-loop pole mass) by a finite shift due to the renormalization procedure. Using the shifted 
mass for the external particle but the tree-level mass for internal particles results in IR di- 
vergences. On the other hand, in general, inserting the shifted chargino mass everywhere 
yields UV divergences. Consequently, we stick to our choice of imposing on-shell conditions 
for the two charginos and one neutralino. The conditions read 



- II 

Ret^oip)] xiip) 
11 



)2— Tn^ - 



These conditions can be rewritten in terms oi six equt 
field renormalization constants or three complex ones 



Re 

Re m 
Re 
Re 

For the further 



= (z = l,2), (163) 
= . (164) 

of six equations defining six real parameters and 

= , 



.Si? ^^2 



Xi 



m 



,0 (sto(m|o) + Efo(m|o)) + Ef^(m|o) + S|o«(m^o; 



Xi 



11 



= 
= 



--^^^ -^-(m|o) + s|^(m|„; 



J 11 



= 



(165) 
(166) 
(167) 
(168) 



determination of the field renormalization constants we also impose 

(z = l,2) , 



lim 



o 9 
X.- 



lim 



+ m-±)[ReV(p)]^_ 

-Xi [p) 



p2 _ ^2 

X, 



+ m^o)[ReS^o(p)]^^ ^ 



p2 - m|o 

Xi 



X\{p) 











(169) 
(170) 
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This leads to the following set of equations: 



Re 



i S|oKo) + S^(m|o; 



Re 

,2 



X ^ X- 



.^.(Ej;Kj) + Eg(m^,„)) 



+m^o (S|o^'(m|o) + S|o^'(m2o 



Re 



XJ' 



Xi' 



S^oKo)-S«oKo) 



11 



11 








(171) 
(172) 

(173) 
(174) 



where we have used again the shorthand S'(m^) = ^|^2_^2- It should be noted that 



Eq. f ll74p is already fulfilled due to the Majorana nature of the neutralinos. Inserting 
Eqs. f ll55p - fll62p for the renormalized self-energies in Eqs. (I165P - (11681) and solving for 
[5M^-]^^ and [5M-j^o] results in 



Re[5M,o]^^ 
Im [SM^o]^^ 



Re 
Re 

-Re 



X ^ X- 



X ^ X. 



X ^ X? 



X ^ X? 



m-± Im 


5z^_ - 5zf: 




J n 2 ^» 




n 



E-(™|.)-E-(m|.: 



Xi \^"x"^"''xi^ ' "x"v"xi- 

SRl^2 \ \^5L/ 2 



J 11 



E|.«(,„4;) - El."!™!; 



11 



^oIm[5Z^o] 



11 



(175) 
(176) 

(177) 
(178) 



where we have used already Eqs. f ll72p and f ll74p . Using Eqs. (115 ip - f ll54p . these conditions 
lead to EOlES] 



5Mi 



2N*^ [N;^S{Mzs^ cos /3) - N*^6{Mzs^ sin , 



A^i*2 [2N*^6iMzc^ cos /3) - 2N*^6iMzc„ sin /3) + iV^s^Ma] 



+ 2iV*3iV*45/i + m^o ( ReE|o (mio) - i Im dZ^o ) + ReS|ct(m'-o 



9 fr/* r/* T/* V* — r/* r/* 1/* V* ^ 

^ l'-^ll'-^22 ^11 ^22 '-^12'-^21 ^12 *^21> 



x"\---xi'' 

X 



■^5L/_2 



11 



(179) 



TT* y* 
^22 ^22 



TT* Y* 



,_2 ^ , T>.^R^^^2_^ 
^~-^L t^2 \ , ^„s:^R 1^2 



m-± (^ReS|4m2±) + ReS^_(m^±) - i Im{(5Z^_ - 5Z^_} ) + 2ReS?f^(m^± 



X ^ xt 

•S'L/,^2 



11 



5/i 



m,. (^ReS|_(m^,) + ReSf.(m^,) - z Im{5Z|_ - 5Z«_}j + 2ReSf (m^,; 
+ 2(f/*2f/2*i - t/nf/22) V^i2^22 5(v^Mh/ Sin /3) 
+ 2f/*2f/22(^i2V^;i - ^1*1^22) K^Mw cos/3) 

1 



J 22 



(180) 



2 (f/rit/|2V;*i^2 - f/l*2f/2*l ^12^2*1^ 
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J J* V* 



-zIm{5Z^_-5Zf_} 
2(f/*2f/2*i - f/nf/22)^i*i^;i 5(v^Mh/ cos (3) 
2f/rit/2i(V^i2^2i - Vl*i^*2)'5(v^MH/sin/3) 



2ReSp( 

A 



22 



ReSj_(m2±) + ReSf_(m2±; 



2ReSf^(m2±^ 

X Xi' 



11 



(181) 

Equations fllTip - ( I173p define the real part of the diagonal field renormahzation constants 
of the chargino fields and of the lightest neutralino field. We generalize the latter result for 
the diagonal field renormahzation constants of the other neutralino fields imposing Eq. fll73p 
also for the components /c = 2, 3, 4 - though we do not define them fully on-shell; see below. 
The imaginary parts of the diagonal field renormahzation constants are still undefined. For 
the definition of the imaginary parts, we use Eqs. (11761) and (I178p . where the latter one is 
generalized for the components fc = 2, 3, 4 and is imposed to hold also for those neutralinos. 
Now, for the charginos and the lightest neutralino Eqs. (11761) and (I178p already define the 
imaginary parts of [(5M^±]^^ [i = 1,2), and [5M^o] which means that a further condition 
for the imaginary part of the field renormahzation constants of the charginos and the lightest 
neutralino is required, and we just set them to zero (see below Eqs. (I185P and (I188P ) which 
is possible as all divergences are absorbed by other counterterms. The off-diagonal field 
renormahzation constants are fixed by the condition that 



ReS^- (p) Xi (p) 



9 2 



ReE^o(p) x'iip) 

J kl 



(^,J = 1,2) 



(fc,/ = 1,2,3,4). 



:i82) 



:i83) 



Finally, this yields for the field renormahzation constants 



Re 


5Z"/« 

X 


a 


Im 


X 


a 




SZ^J^ 






X 




Re 


'sz-o'_ 


kk 



Im 



5Z,o 



5Z,o 



kk 



kl 



-Re 



E^/^(m2±) 

X Xi 



-Re<^S 



SL 



^ Xi' 



lm6M, 



2 2 

ml± — ml± 

Xi Xj 



Re 



Xf X ^ Xf^ Xi Xj X ^ X 



+ m-±S-_ (m- 
-Re 



+ m-±J:l^^^^(mi±) - m-± 6M^^ - m-± 6M 



X 



Xi/j 



E^o(m|o; 



(184) 

- n 

(185) 
(186) 



^3 

(187) 



Xk 



Xk 

i 



Sfo(m2 



Xl' 



Xk 



E|fK.)+E?,r(™|.: 



kk 



Re {J:f^(mio) - ^%l^(mio 



SL, 



Xk 



Xk' 



Im (5M^o 



k=l 



0, 



(188) 



J kk 



Xk 



Xl 



Re 



xr 



xr' 



Xl 
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+ m^oS|,f (m|o) + m^oS|,?^(m|o) - m^o6M^o - m^o^Mto 



J kl 



The Eqs. f ll63p . fll64p result in three on-shell masses in the neutralino/chargino sector. 
Therefore the three neutralino masses 1^X234^ '^^ other hand, require a finite shift for 
their on-shell value. We have checked that this shift (for the scenarios under investigation 
in Sec. H]) is numerically small and does not change our results. Consequently, these shifts, 
though formally necessary, are not further taken into account to simplify the numerical 
evaluation. 



The field renormalization constants for squark, quark, gluino, gauge boson as well as 
chargino and neutralino fields that have been derived in Sec. 12.11 - 12.31 are constructed via 
the multiplicative renormalization procedure in a symmetry conserving way, absorb the diver- 
gences accordingly and are defined via on-shell renormalization conditions. In the presence 
of complex phases and nonvanishing absorptive parts of the self-energy type corrections, 
further wave function corrections may arise that are not part of the field renormalization 
constant but can be taken into account by additional Z factors; see the Appendix. 



3 Calculation of loop diagrams 

In this section we give some details about the calculation of the higher-order corrections 
to the partial decay widths of scalar quarks. Sample diagrams are shown in Figs. [U - El 
Not shown are the diagrams for real (hard and soft) photon and gluon radiation. They are 
obtained from the corresponding tree-level diagrams by attaching a photon (gluon) to the 
electrically (color) charged particles. The internal generically depicted particles in Figs. [T]- 
[7] are labeled as follows: F can be a SM fermion /, chargino x^, neutralino xt^ gluino g; 
S can be a sfermion /j or a Higgs boson hn] V can be a photon 7, gluon g, or a massive SM 
gauge boson, Z or W^. For internally appearing Higgs bosons no higher-order corrections 
to their masses or couplings are taken into account; these corrections would correspond to 
effects beyond one-loop order@ For external Higgs bosons, as described in Sec. 12. 2[ the 
appropriate Z factors are applied and on-shell masses (including higher-order corrections) 
are used. 

Also not shown are the diagrams with a gauge/ Goldstone-Higgs boson self-energy contri- 
bution on the external Higgs boson leg. They appear in the decay ^2 iihn, Fig. [H with a 
Z/G-hn transition and in the decay ^2 hH^, Fig. El with a /G^-H'^ transition^ The 
corresponding self-energy diagram belonging to the process t2 tiZ or ^2 — ^ h^^, respec- 
tively, yields a vanishing contribution for external on-shell gauge bosons due to e ■ p = for 

= M| (p^ = M^), where p denotes the external momentum and e the polarization vector 
of the gauge boson. 

Furthermore, in general, in Figs. [T]-[7]we have omitted diagrams with self-energy type 
corrections of external (on-shell) particles. While the contributions from the real parts of the 
loop functions are taken into account via the renormalization constants defined by on-shell 

We found that using loop corrected Higgs boson masses in the loops leads to a UV divergent result. 
19 pj-om a technical point of view, the /G'^-H'^ transitions have been absorbed into the respective 
counterterms, while the Z /G-hn transitions have been calculated explicitly. 
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Figure 1: Generic Feynman diagrams for the decay t2 — ^ ^i^^n (^^ = 1, 2, 3). F can be a SM 
fermion, chargino, neutralino, or gluino; 5" can be a sfermion or a Higgs boson; V can be a 
7, Z, W^, or g. Not shown are the diagrams with a Z-hn or G-hn transition contribution 
on the external Higgs boson leg. 



renormalization conditions, the contributions coming from the imaginary part of the loop 
functions can result in an additional (real) correction if multiplied by complex parameters 
(such as At). In the analytical and numerical evaluation, these diagrams have been taken into 
account via the prescription described in the Appendix. The impact of these contributions 
will be discussed in Sec. IH 

Within our one-loop calculation we neglect finite width effects that can help to cure 
threshold singularities. Consequently, in the close vicinity of those thresholds our calculation 
does not give a reasonable result. Switching to a complex mass scheme [OS] would be another 
possibility to cure this problem, but its application is beyond the scope of our paper. 

Finally it should be noted that the purely loop induced decay channels t2 ill / 9 yield 
exactly zero due to the fact that the decay width is proportional to e-p and the photon/gluon 
is on-shell, i.e. e ■ p = 0. 

The diagrams and corresponding amplitudes have been obtained with FeynArts p7] . 
The model file, including the MSSM counterterms, is largely based on Ref. [50], however 
adjusted to match exactly the renormalization prescription described in Sec. [2l The further 
evaluation has been performed with FormCalc [68] . 
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Figure 3: Generic Feynman diagrams for the decay t2 tg. F can be a SM fermion, 
chargino, neutralino, or gluino; S can be a sfermion or a Higgs boson; V can be a 7, Z, W^, 
or g. 
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Figure 4: Generic Feynman diagrams for the decay ^2 — ^ tx^ {k = 1,2,3,4). F can be a 
SM fermion, chargino, neutralino, or gluino; S can be a sfermion or a Higgs boson; V can 
be a 7, Z, VT^*", or g. 




Xj Xj Xj Xj 



b 




Figure 5: Generic Feynman diagrams for the decay t2 — )■ bxf {j = 1,2). F can be a SM 
fermion, chargino, neutrahno, or gluino; S can be a sfermion or a Higgs boson; V can be a 
7, Z, or 5f. 



Ultraviolet divergences 

As regularization scheme for the UV divergences we have used constrained differential renor- 
malization [69], which has been shown to be equivalent to dimensional reduction [70] at the 
one-loop level [68]. Thus the employed regularization scheme preserves SUSY [7T],[72] and 
guarantees that the SUSY relations are kept intact, e.g. that the gauge couplings of the 
SM vertices and the Yukawa couplings of the corresponding SUSY vertices also coincide to 
one-loop order in the SUSY limit. Therefore no additional shifts, which might occur when 
using a different regularization scheme, arise. All UV divergences cancel in the final result. 
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Figure 6: Generic Feynman diagrams for the decay ^2 hiH^ {i = 1,2). F can be a SM 
fermion, chargino, neutralino, or gluino; S can be a sfermion or a Higgs boson; V can be a 7, 
Z, W^, or g. Not shown are the diagrams with a W~^-H^ or G'^-H^ transition contribution 
on the external Higgs boson leg. 
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Figure 7: Generic Feynman diagrams for the decay ^2 hiW^ (i = 1,2). F can be a SM 
fermion, chargino, neutralino, or gluino; S can be a sfermion or a Higgs boson; can be a 
7, Z, or g. 
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Infrared divergences 



The IR divergences from diagrams witli an internal plioton or gluon have to cancel with 
the ones from the corresponding real soft radiation. In the case of QED we have included 
the soft photon contribution following the description given in Ref. [61]. In the case of 
QCD we have modified this prescription by replacing the product of electric charges by the 
appropriate combination of color charges (linear combination of Ca and Cp times a^). The 
IR divergences arising from the diagrams involving a 7 (or a g) are regularized by introducing 
a photon (or gluon) mass parameter, A. While for the QED part this procedure always works, 
in the QCD part due to its non-Abelian character this method can fail. However, since no 
triple or quartic gluon vertices appear, A can indeed be used as a regulator (the appearance 
of the non-Abelian gluino-gluino-gluon vertex does not pose a problem here [12] )• All IR 
divergences, i.e. all divergences in the limit A — >■ 0, cancel once virtual and real diagrams for 
one decay channel are added. 

Special care has to be taken in the decay modes involving scalar bottom quarks. Using 
tree-level sbottom masses yields a cancellation of IR divergences to all orders for all ^2 
decay modes. However, inserting the one- loop corrected sbottom masses (see Sec. I2.1.ip . as 
required for consistency, we found cancellation to all orders of the related IR divergences, 
except for the decay modes ^2 ~^ biW~^. Within these decays the tree- level relation required 
by the SU{2) symmetry Mq^{t) = Mq^{h), corresponding to 



I Oil I bi 



I £>12 I 62 



\UlJrnl + \Ui^^ 



2 2 _2 , _2 ,^2 
12 



+ m;- cos 2/3 



(190) 



has to be fulfilled to yield a cancellation of all IR divergenceso On the other hand, the 
requirement of on-shell sbottom masses as well as an intact SU (2) relation at the one-loop 
level leads to the necessity of a shift in the scalar bottom masses; see Eq. f l48p . Therefore 
Eq. fll90p is "violated" at the one- loop level, introducing a two- loop IR divergence in r(t2 — ^ 
biW~^). In order to eliminate this two-loop IR divergence we introduced a counterterm in 
the t2biW vertex. 



If/ |2^2 1^ ,2 2 
I Oil I fei ' 012 1 b2 



I oil I bi ' " 
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^12 1 62 



shift 



X IR div 



(191) 



to restore the tree-level SU{2) relation. The left term in Eq. fll9ip contains only "tree- 
level" values, while the index "shift" refers to inserting the one-loop masses and mixing 
matrices. The IR divergence has been taken from Eq. (B.5) of Ref. [73] (it can also be found 
in Ref. [H]), and reads (in our case) 



IR div 



a 



2 xt ln{xt)_ f mi^Mw 



"^t2(l 



A2 



Xb In(xb) ^^(miMw 



A2 



(192) 



with 



^1 - 4m,-^Miy/(m? + ^0 - {Mw - mi^f) - 1 
y/l -Ami^Mw/{ml^ + ^0 - {Mw - m,-J2) + l 
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(193) 



Equation has been deduced via M?^ = It/qupm?^ + |f7giJ2m| - Mlc2fs{ll - Q^sl) 
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/I -Am~, Mw/{ml + zO - [Mw - m-,f) - 1 

Xb = \ , (194) 

^1 -Am-,Mw/iml + tO - (M^ - m~,f) + 1 

where denotes an infinitesimally small imaginary part. After including this tree- level 
relation restoring counterterm we find an IR finite result to all orders as required. 

We have furthermore checked that our result does not depend on AE defining the energy 
cut that separates the soft from the hard radiation. Our numerical results have been obtained 
for = 10"^ X mz. 



Tree-level formulas 



For completeness we show here also the formulas that have been used to calculate the tree- 
level decay widths: 



\C{h,h,h^)\n'/\ml,ml,ml 



16 vrmf 

12 



(n = l,2,3) 



\C{h,ii,Zr\^/\mlmlM, 



16 TT Ml m? 
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(A; = 1,2,3,4) 



r*-^^~^ ^ bxt) = [ {\Cih. xt)L\' + |C7(t2, 6, x;)/?!') (m| - - m^^) 
- 4 Re{C(t2, b, xt)! C{h, b, x+)k} m- 



16 vrmf 
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(z = l,2) 
(^ = 1,2) : 
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where X{x,y,z) = {x — y — z)^ — Ayz and the couphngs C(a,6, c) can be found in the 
FeynArts model files [73]. C{a,b,c)L,R denote the part of the coupling that is proportional 
to(]lT75)/2. 



4 Numerical analysis 

In this section we present a numerical analysis of all 15 decay channels. In the various figures 
below we show the partial decay widths and their relative correction at the tree-level ( "tree" ) 
and at the one- loop level ("full"), 

(202) 
(203) 

(204) 

where xy denotes the specific final state. The total decay width is defined as the sum of all 
15 partial decay widths, 

"pfuU "ptree 

rr = E xy) , T'^l ^ Yl r'""(^~2 xy) , 5r,o,/rJr = ^^W^ • (205) 

We also show the absolute and relative changes of the branching ratios. 



■ptree 


= r*'-^^(t2 ^ xy) 


■pfuU 


^ T'-\h ^ xy) 


^p/ptree 


pfuU ptrcc 


ptree ' 



^ ^"^^L?"^^ ^ (206) 



tot 



BR'-'' = ^ , (207) 



r^""(t2 ^ xy) 

pfuU 



tot 



Topfull _ TOO tree 

<5BR/BR = . (208) 

The last quantity is crucial to analyze the impact of the one-loop corrections on the phe- 
nomenology at the LHC and the ILC. 

4.1 Parameter settings 

The renormalization scale hr has been set to the mass of the decaying particle, i.e. fin = m^^. 
The SM parameters are chosen as follows; see also 

• Fermion masses (on-shell masses, if not indicated differently) : 

me = 0.51099891 MeV , m^^ = MeV , 

= 105.658367 MeV , = MeV , 

rrir = 1776.84 MeV , m^, = MeV , 



Using the most up-to-date values from Ref. [75] would have a negligible impact on our numerical results. 
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m„ = 53.8 MeV , = 53.8 MeV , 

mc = 1.27 GeV , m, = 104 MeV , 

mt = 171.2 GeV , mf,(mb) = 4.2 GeV . (209) 

According to Ref. [H], m^, is an estimate of a so-called "current quark mass" in the 
MS scheme at the scale /i ~ 2 GeV. and rrih are the "running" masses in the MS 
scheme. The top quark mass as well as the lepton masses are defined OS. m„ and 
are effective parameters, calculated through the hadronic contributions to 

f=u,c,d,s,b J 

• The CKM matrix has been set to unity. 

• Gauge boson masses: 

Mz = 91.1876 GeV , Mw = 80.398 GeV . (211) 

• Coupling constants: 

a = — = 1/137.0359895 , as{Mz) = 0.1176 , (212) 
47r 

where the running and decoupling of can be found in Sec. 12.1.31 

The Higgs sector quantities (masses, mixings, etc.) have been evaluated using FeynHiggs 
(version 2.6.5) [MH3Z]@ 

We will show the results for some representative numerical examples. The parameters 
are chosen according to the two scenarios, SI and S2, shown in Tab. [H but with one of the 
parameters varied. The scenarios are defined such that all decay modes are open simulta- 
neously to permit an analysis of all channels, i.e. not picking specific parameters for each 
decay. We will start with a variation of m^^, and show later the results for varying (p^^. 
The scenarios are in agreement with the MSSM Higgs boson searches at LEP [771[78]. Too 
small values of the lightest Higgs boson mass would be reached for tan /3 ^ 9.4 (4.6) within 
SI (S2) as given in Tab. [T][^ In order to avoid completely unrealistic spectra, the following 
exclusion limits HH hold in our two scenarioO: 



rrii^ > 95 GeV, m^^ > 89 GeV, nig > 379 GeV, ms, > 73 GeV, 

m^o > 46 GeV, m^± > 94 GeV, nig > 308 GeV. (213) 



As default value within FeynHiggs, jiR — rrit is used. Furthermore we have neglected the (in our case 
small) corrections of 0{abas,atab.,a^) (via a small modification in the code) and used the top pole mass 
for the evaluation of the MSSM Higgs boson sector quantities. 

^■^ While in these scenarios we are not aiming to yield a light Higgs boson mass value around ~ 125 GeV, 
it should be noted that such a value is in principle in agreement with not too heavy scalar top quarks [79j 
that can be produced via e+e~ — > t\i2 at the ILC(IOOO). 

The relatively light scalar quark masses and especially the light gluino mass are potentially in conflict 
with recent SUSY searches at the LHC [8^ (although no fully model-independent results have been pub- 
lished). However, as stressed above, the parameters are chosen to be able to analyze as many decay modes as 
possible simultaneously. For a realistic collider analysis these bounds |80| will have to be taken into account. 
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Seen. 


tan /3 


Mh± 






^1, 




At 


^6 


Ml 


M2 


M3 


SI 


20 


150 


650 


0.4 


0.7 m-t^ 


200 


800 


400 


200 


300 


350 


S2 


20 


180 


1200 


0.6 mi^ 


0.8 m-t^ 


300 


1800 


1600 


150 


200 


400 



Table 1: MSSM parameters for the initial numerical investigation; all parameters (except 
of tan/3) are in GeV. We always set m^^(mfe) = 4.2 GeV. In our analysis Mq^{= Mj^), 
Mi^{= Muj^ = Mcj^), and Mi^{= Mj^ = M^^ = Mj^) are chosen such that the values of m^-^, 
771^2 , and m^_^ are realized. For the b sector the shifts in Mq^^{h) as defined in Eqs. (148|) and 
(ISO]) are taken into account. The values for At and Af, (= A,-) are chosen such that charge- 
or color-breaking minima are avoided [ST] . 



Seen. 


tan/3 










SI 


2 


260.000 


650.000 


305.436 


455.000 


20 


260.000 


650.000 


333.572 


455.000 


50 


260.000 


650.000 


329.755 


455.000 


S2 


2 


720.000 


1200.000 


769.801 


960.000 


20 


720.000 


1200.000 


783.300 


960.000 


50 


720.000 


1200.000 


783.094 


960.000 



Table 2: The stop and sbottom masses in SI and S2 and at different tan (3 for the numerical 
investigation; all masses are in GeV and rounded to 1 MeV. 



A few examples of the scalar top and bottom quark masses in SI and S2 are shown in 
Tab. [2j The values of m^^ allow copious production of the heavier stop at the LHC. For other 
choices of the gluino mass, > mj^, which would leave no visible effect for most of the 
decay modes of the ^2, the heavier stop could also be produced in gluino decays at the LHC. 
Furthermore, in SI (even for the nominal value of m^^ as given in Tab. H]) the production 
of ^2 at the ILC(IOOO), i.e. with i/i = 1000 GeV, via e+e" i\i2 will be possible, with all 
the subsequent decay modes (III) - (I7j) being open. The clean environment of the ILC would 
permit a detailed study of the scalar top decays. For the lowest values shown in the plots 
below, mi^ > 570 GeV, we find (via a tree-level calculation) cr(e+e~ — )■ i\t2) ~ 1.5 fb, i.e. an 
integrated luminosity of ~ 1 ab~^ would yield about 1500 ^2- This number drops to ~ 280 ^2 
for the masses shown in Tab. [2 The ILC environment would result in an accuracy of the 
relative branching ratio (Eq. f l208p ) close to the statistical uncertainty: a BR of 30% could 
be determined to ~ 5% for the lowest m^^ values and to about 11% for the values given 
in Tab. [2l Depending on the combination of allowed decay channels a determination of the 
branching ratios at the few percent level might be achievable in the high-luminosity running 
of the ILC(IOOO). 
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The numerical results we will show in the next subsections are of course dependent on 
choice of the SUSY parameters. Nevertheless, they give an idea of the relevance of the full 
one-loop corrections. As an example, the largest decay width is T{t2 — )■ tg), dominating the 
total decay width, Ftot, and thus the various branching ratios. For other choices of rrig with 
rrig > the corrections to the decay widths would stay the same, but the branching ratios 
would look very different. Channels (and their respective one-loop corrections) that may 
look unobservable due to the smallness of their BR in the plots shown below, could become 
important if other channels are kinematically forbidden. 

4.2 Full one-loop results for varying m^^ 

The results shown in this and the following subsections consist of "tree", which denotes 
the tree-level value and of "full", which is the partial decay width including all one-loop 
corrections as described in Sec. [31 We start the numerical analysis with partial decay widths 
of t2 evaluated as a function of m^^, starting at m^^ = 570 GeV up to m^-^ = 3 TeV, which 
roughly coincides with the reach of the LHC for high- luminosity running. The upper panels 
contain the results for the absolute value of the various partial decay widths, r(t2 xy) 
(left) and the relative correction from the full one-loop contributions (right). The lower 
panels show the same results for BR(t2 — ^ xy). 

Since in this section all parameters are chosen to be real, no contributions from absorptive 
parts of self-energy type corrections on external legs can contribute. This will be different 
in Sec. |01 

In Figs. |H]-[ID]we show the results for the process t2 — )■ tihn {n = 1, 2, 3) as a function of 
■ These are of particular interest for LHC analyses [8l[9] (as emphasized in the Introduc- 
tion). The dips at m^-^ ~ 819, 948, 971, 1264, 1303 GeV (for all three figures) in the scenario 
SI are effects due to the thresholds rrit + m^o ^ 3 ^ = "^ii and nit + f^g = fTiii (in this order) of 
the self-energies Sj^^^il'^f ) in the renormalization constants [^Z^],^^^!) ^^t) and 5m? . One 
can see that the size of the corrections of the partial decay widths is especially large very 
close to the production thresholcQ from which on the considered decay mode is kinematically 
possible. Away from this threshold relative corrections of ~ -|-10%, —20%, —5% are found 
for hi, h2, h^, respectively. In (all) the plots the value of m^^ which m^^ +"^42 ~ 1000 GeV 
is shown as a vertical line, i.e. the region where the heavier stop can be produced at the 
ILC(IOOO). In these regions the size of the corrections amounts up to ~ +20%, —10%, +10% 
for the three neutral Higgs bosons. The BRs are at the few percent level for all three channels 
for the two numerical scenarios. The relative change in the BRs for the masses accessible at 
the ILC(IOOO) are about +8%, —21%, —1% for hi, /i2, h^, respectively. For lager masses, 
only accessible at the LHC, the one-loop corrections are around +10%, —25%, and —5%. 
Depending on the MSSM parameters (and the channels kinematically allowed) the one-loop 
contributions presented here can be relevant for analyses at the ILC as well as at the LHC. 

Next, in Fig. [TT]we show results for the decay T(t2 — >■ iiZ). The dips due to the thresholds 
in [(5Zf]_^^2i) ^^t! and Sm"^^ are the same as before. The relative corrections to the partial 

It should be noted that a calculation very close to the production threshold requires the inclusion of 
additional (nonrelativistic) contributions, which is beyond the scope of this paper. Consequently, very close 
to the production threshold our calculation (at the tree- or loop-level) does not provide a very accurate 
description of the decay width. 
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decay width in SI range between +8% at low m^^, i.e. in the "ILC(IOOO) regime", to —5% 
at large m^^; with the exception of the region close to thresholds. Within S2 the relative 
corrections stay below ~ 5%. The BR(t2 — ^ tiZ) is larger than 15% for the smallest m^^ 
values in the two scenarios. This drops below 3% for m^^ ^ 2.5 TeV. The relative change 
for masses accessible at the ILC(IOOO) is found at the few percent level. 

The results for the decay t2 — )■ tg are presented in Fig. [121 We see that for the relative 
corrections of the partial decay width up to 48% (22%) are reached for m^^ = 570 (980) GeV 
in SI (S2), i.e. at the smallest possible value and decrease for increasing m^-^. It should 
be noted that in this case the hard and soft QCD radiation can be very large and the two 
compensate each other. The BR turns out to be very large and growing with m^-^, where 
values larger than 50% are found. Within SI the relative corrections can reach up to +20% in 
the production threshold region and are larger than +12% in the parameter space accessible 
at the ILC(IOOO). For large m^-^ the corrections range between +11% and +15% in the two 
scenarios. 

Now we turn to the decays t2 txt {k = 1,2,3,4), with the results shown in Figs. [13] 
^ [T6l Since fi, Mi, and M2 are roughly of the same order, the four states are a mixture of 
gauginos and Higgsinos. Consequently, the partial decay widths are found to be roughly the 
same. The larger partial decay widths for the decay modes ^2 — ^ ^Xfc with k = 1,2 {k = 3,4) 
are found in SI (S2) and are ~ 15 — 25 GeV. For SI we find relative one-loop corrections 
ranging between 0% and ~ —30% for T{i2 — )■ tXi) — ^ 'tX2)y where the smaller values 
are reached for small m^-^. T{t2 — > tXs) receives one-loop contributions between +8% and 
—16%, while for T{i2 — > tx^) we find —18% to —35% with the exception of very small m^^, 
where the partial decay width itself is negligible. Within S2 the corrections stay at the few 
percent level for 1(^2 tx?), while for 1(^2 ^ tx"), ^(h ^ ^Xs) and 1(^2 ^ txl) they 
range between —10% and —30%. Following the size of the partial decay widths, also the 
branching ratios are roughly the same for the four decay modes. The relative changes in 

the BRs for m^-^ + m^^ < 1000 GeV are 15%, -18%, -6%, and -30% for k = 1,2,3,4, 

respectively. Especially, for this parameter range, for the decays to the two lighter neutralinos 
the branching ratios are ~ 5% and ~ 10%, i.e. the one-loop corrections can be crucial to 
match the anticipated ILC precision. 

Next in Figs. [I7l[l8]we present the results for ^2 bxf (j = 1,2). The size of the partial 
decay widths and branching ratios for ^2 bxf (^2 — ^ ^X2^) are roughly the same as for 
^2 txl with k = 1,2 {k = 3,4). For r(t2 — > bxt) we find relative corrections starting at 
—5% at low 771^2 in SI down to ~ —35% at high m^-^ in both scenarios. The partial decay 
width r(t2 — ^ bx2) is very small in SI for m^^ + m^^ < 1000 GeV. Because of this smallness, 
and additionally pronounced due to the vicinity of the production threshold, the relative size 
of the corrections becomes huge and is not reliable anymore. For higher m^-^ values we find 
relative corrections between —20% (—10%) to —30% in SI (S2). A large branching ratio of 
~ 14% in the ILC(IOOO) accessible regime is reached in SI in the decay t2 — bxt , where the 
one-loop corrections are ~ —20%. Again the one-loop corrections can be crucial to match 
the ILC precision. 

We now turn to the decay modes ^2 ~^ biH^ {i = 1,2). Results are shown in Figs. [T9l 
[2U| (which have been used for the investigations in Ref. [SSj). In Fig. [T^ several peaks and 
dips can be observed. Within SI the first (fourth) peak at m^-^ ~ 571 (638) GeV stems 
from the threshold m^-^ + Mh± (Mw) = in the self-energies ^^^^^^^(m? ) entering the 
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renormalization constants [^Z^j^^^^i- The second dip at m^^ ~ 596 GeV comes from the 
threshold rrig + = m^^ . The third and the fifth dip come from the same thresholcQ rrit + 
m-± = m^^ at m^^ ~ 601, 823 GeV, respectively^ The sixth dip at m^^ ~ 1282 GeV comes 
from the threshold rrit + m-± = m^^. Within S2 the peak/dip at m^^ ~ 1130 GeV is the 
threshold m^^ + Mw = fn^ . The other peaks/dips do not appear as the values for the stop 
masses are different. Also in Fig. [20] some peaks/dips appear. Within SI the first peak/dip is 
visible at m^^ ^ 571 GeV, due to the threshold mi^+MH± = rn^^ (in the self-energy Sjj^^(m? ) 

entering the renormalization constant [^Z^j^-j^). Within S2 the "apparently single" dip is in 
reality two dips at m^^ ^ 1163 (1164) GeV coming from the thresholds m^^ {rriH) = rrij^^ 
(it should be noted that the internal Higgs boson masses are tree- level masses). 

The absolute value of the partial decay widths is relatively small, staying below ~ 
1.2 (0.2) GeV for 1(^2 hH+) (1(^2 ^ W^)). The relative size of the one-loop cor- 
rections to r(t2 biH+) ranges between 12% and -27% (-20%) in SI (S2). For 

r{t2 — )■ &2-f^^) very large corrections are found for the smallest m^^ values, dropping to 
values close to zero for larger masses. Because of the small partial decay widths also the 
branching ratios are at or below the 1% level. Only if other channels were kinematically sup- 
pressed, these decays could play a relevant role, and the one-loop effects could be expected 
at the level of one-loop contributions to the partial decay widths itself. 

Results for the other decay modes involving scalar top and bottom quarks, ^2 biW~^ 
{i = 1,2), are shown in Figs. [211 Ell Also these decay modes have been analyzed in detail 
in Ref. [3H]. The peaks/dips are the same ones as for the decays t2 — )■ hiH^ . On top of that 
due to different renormalization constants entering the calculation one observes the following: 
within SI peaks appear at m^^ ~ 618 (656,657) GeV due to m^^ -|- Mz {mA^rnu) 
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and at m^-^ ^721 GeV due to rrit + = ^t2- The "knee" at m^-^ ~ 1303 GeV results 
from rrit + = i^ii ^?i2('^fi) entering the renormalization constant 6Yt. Within S2 
one hardly visible dip can be found at ~ 1039 GeV from m^^^ -|- Mz = ^i^- The 
"knee" at mi^ ~ 1163 (1164) GeV is the same one as in Fig. [20l The absolute size of 
r(t2 — ^ hiW^) is found to be between ~ 5 GeV and ~ 25 GeV, depending on m^^ and the 
scenario. r(t2 i>2^^)i on the other hand, is found to be tiny for nearly all m^-^ values. 
However, the smallness of r(t2 h2W^) is a purely coincidental effect. A slightly smaller 
would yield a width of (9(0.1 GeV). The relative corrections to r(t2 hiW^) range 
between +10% (+15%) and -5% (0%) for SI (S2). For r(t2 ^ ^2^^+) we find corrections of 

10% (-30%) to +18% (+3%) in SI (S2), where it has to be kept in mind that the partial 

decay width itself is tiny in our scenarios SI and S2. The branching ratio for t2 — hiW^ can 
reach up to ~ 18% in the parameter range with m^^ + ^ 1000 GeV. Here the relative 
one-loop effect on the BR is ~ —5% and could be important to reach the ILC precision. 



It should be remembered that mg^ changes its value when the value of m^^ is changed. 

For these two different input parameters (m^.^ w 601, 823 GeV) we get coincidcntally m^^ sa 349.14 GeV. 
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Figure 8: T{t2 — )■ tihi). Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^-^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 9: T{t2 — )■ ^1/^2)- Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^-^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 10: T{t2 — )■ tih^). Tree- level and full one- loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 11: r{t2 ^ tiZ). Tree- level and full one- loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^_^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 



42 






1 1 


1 1 










/ X 







/ / ' / 


SI, tree - 






SI, full 






S2, tree " 




1 1 


S2, full 



5BR/BR 
22% 



500 



1000 



1500 



2000 



2500 



3000 




3000 



Figure 12: T{t2 — )■ tg). Tree- level and full one- loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^_^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 



43 





Figure 13: r(t2 — ^ tXi)- Tree- level and full one- loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^-^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 14: T{t2 — )■ tx2)- Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^-^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 15: r(t2 — ^ tXs)- Tree- level and full one- loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^-^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 16: r(t2 — ^ tXi)- Tree- level and full one- loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^-^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 17: T{t2 — )■ bxt)- Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 18: T{t2 — > bxt)- Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^_^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 19: T{t2 — )■ hiH^). Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 20: T{t2 — )■ h2H^). Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 21: T{t2 — ?■ biW~^). Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^_^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 
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Figure 22: T{t2 — ?■ b2W~^). Tree-level and full one-loop corrected partial decay widths are 
shown with the parameters chosen according to SI and S2 (see Tab. [1]), with m^-^ varied. The 
upper left plot shows the partial decay width; the upper right plot shows the corresponding 
relative size of the corrections. The lower left plot shows the BR; the lower right plot shows 
the relative correction of the BR. The vertical lines indicate where m^-^ + m^_^ = 1000 GeV, 
i.e. the maximum reach of the ILC(IOOO). 



53 



4.3 Full one-loop results for varying (f^^ 

In this subsection we analyze the various partial decay width@ and branching ratios as a 
function of (pAf The other parameters are chosen according to Tab. [1] Thus, within SI 
we have m^-^ + mf^ = 910 GeV, i.e. the production channel e~^e~ — )■ i\i2 is open at the 
ILC(IOOO). Consequently, the accuracy of the prediction of the various partial decay widths 
and branching ratios should be at the same level (or better) as the anticipated ILC precision. 
It should be noted that the tree-level prediction depends on ip^^ via the stop mixing matrix. 

When performing an analysis involving complex parameters it should be noted that the 
results for physical observables are affected only by certain combinations of the complex 
phases of the parameters fi, the trilinear couplings At, Ab, . . . , and the gaugino mass param- 
eters Ml, M2, M3 [S311H2]- It is possible, for instance, to rotate the phase ipM2 away. Exper- 
imental constraints on the (combinations of) complex phases arise, in particular, from their 
contributions to electric dipole moments of the electron and the neutron (see Refs. [83|[8^ 
and references therein), of the deuteron [85] and of heavy quarks [86]. While SM contribu- 
tions enter only at the three-loop level, due to its complex phases the MSSM can contribute 
already at one-loop order. Large phases in the first two generations of sfermions can only 
be accommodated if these generations are assumed to be very heavy [ST] or large cancel- 
lations occur [SB]; see, however, the discussion in Ref. [52] • A recent review can be found 
in Ref. [90]. Accordingly (using the convention that (pM2 = 0, as done in this paper), in 
particular, the phase is tightly constrained [91], while the bounds on the phases of the 
third generation trilinear couplings are much weaker. The phases of fi, At, and Ai, enter 
only in the combinations {(fiAt t + V^m) different combinations together with phases of 

Ml or M3). Setting ip^ = (see above) as well as ipMi = fg = ^ (we do not consider these 
phases in this paper) leaves us with ipAt and ipA^ as the only complex valued parameters. 
The dependence on cpA^ on the partial decay widths involving scalar bottom quarks has been 
analyzed in detail in Ref. [38]. Consequently, we focus on a complex At and keep Af, real. 

Since now a complex At can appear in the couplings, contributions from absorptive parts 
of self-energy type corrections on external legs can arise, and their impact will be discussecj^. 
The corresponding formulas for an inclusion of these absorptive contributions via finite wave 
function correction factors can be found in the Appendix. 

As before we start with the decays to Higgs bosons, ^2 — ^ iihn {n = 1,2,3) shown in 
Figs. [23] - [251 The arrangement of the panels is the same as in the previous subsection. 
In Fig. [231 where the partial decay width T(t2 — iihi) is given as a function of ipAt^ oiie 
can see that the size of the corrections to the partial decay width vary substantially with 
ifAf The one-loop effects range from +21 (+16)% to +6 (+4)% in SI (S2). The effect 
of the absorptive parts of self-energy type corrections on external legs (called "absorptive 
contributions" from now on) are at the few percent level. For ipAt = 0, n, 2ti these effects 
vanish (by construction). 

It should be kept in mind that the parameters are chosen such that e"'"e~ — )■ ^^^2 is 
kinematically possible at the ILC(IOOO) in SI, where the knowledge of such a large variation 
can be very important. Also for ^2 — ^ ^1^2, shown in Fig. (211 the variation with ipAt is very 

Again we note, that we do not investigate the decays of here, which would correspond to an analysis 
of CV asymmetries, which is beyond the scope of this paper. 

In a slight abuse of the language "full" still refers to corrections without absorptive contributions. 
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large, ranging from —6% to —24%, again with a non-negligible shift from the absorptive 
contributions. The wiggles in the size of the relative corrections to the partial decay width 
is a result of small numerical variations and not visible in the upper left panel showing the 
full decay width. These variations are enhanced due to the smallness of the tree-level partial 
width; see Eq. f l204p . The results for t2 — )■ tih^ can be found in Fig. I25l Also here the 
size of the corrections shows a large variation with ip^^ , again with non-negligible absorptive 
contributions. Within S2 for real and negative At the partial width becomes extremely 
small at tree-level, leading to (formally) very large relative one-loop corrections. The one- 
loop effects on the branching ratios also vary strongly with ipA^ , following the same pattern as 
the partial decay widths. Effects up to ±8% are reached for BR(t2 — >■ ^i^i), while the other 
two decay modes reach large corrections only where the BRs are relatively small, ^ 1%. 
The one-loop corrections to T(t2 — )■ ti/ii), however, can easily exceed the ILC precision. 

In Fig. I2S] we present the phase dependence for the decay mode t2 — tiZ. While in SI 
the effect of the one- loop corrections to T{t2 — )■ tiZ) varies from ~ +8% to ~ +5%, within 
S2 only a very small variation can be observed. These numbers change if the absorptive 
contributions are taken into account. In both scenarios substantially larger variations are 
found. Within SI (S2) the branching ratio varies with ip^t between ~ 11% and ~ 15% (9.5% 
and 11%). Again the variation of the relative correction of the BR increases visibly via the 
inclusion of the absorptive contributions. The relative corrections reach —4% (—9%) and 
are relevant to match the ILC precision in SI. 

Next we show the results for t2 — )■ tg in Fig. |271 In both numerical scenarios we find a 
substantial variation of the one-loop effects with ipAf The effects range from +28% (+20%) 
to +36% (+24%) in SI (S2), where the effect of the absorptive contributions remains rel- 
atively small. The branching ratio varies strongly with (fAt in SI, ranging from ~ 38% to 
~ 25%, while in S2 it is larger and varies less, being around ~ 41%. The one-loop corrections 
in SI vary between +14% and +18% and are important for physics at the LHC and the ILC. 
Within S2 they are found to be ~ 8%. 

In Figs. [2S] - ED we present the variation of T{i2 — )■ txl) {k = 1,2,3,4) as a function 
of (fAf As for the variation with m^^ also here for = 1, 2 (fc = 3, 4) larger values of the 
partial decay width are found in SI (S2) with a similar size as before. The one-loop effects 
on r(t2 — > txi) for (fAt = are relatively small, at the +3% (—3%) level in SI (S2). The 
variation of (fAt, however, now yields one-loop corrections up to ~ +10% (—10%) in the 
two scenarios, with a small shift induced by the absorptive contributions. r(t2 tx2) also 
exhibits a strong variation with ipAt, ranging from —12% (—18%) to —6% in SI (S2). The 
absorptive contributions in this case can change the result strongly, leading especially in S2 
to a substantially different shape. For r(t2 — ^ ^Xs) the variation in S2 is small. Within SI, 
however, the effects for ipAt = are at the ~ +5% level, while they reach nearly —15% for 
intermediate ipAf The absorptive contributions are small. The last partial decay width of 
the four decay modes, T{i2 — )■ txf), shows a large variation at the one-loop level of nearly 
—20% in SI, where, however, the partial width itself is very small. For (pAt ~ v" a width of 
~ 2 GeV is reached with a variation at the —10% level. Within S2 T{t2 — tx^) varies around 
~ 7 GeV with a one-loop variation between —12% and —17%. The absorptive contributions 
lead to a result smaller by a few percent. Within SI, i.e. with the ILC(IOOO) accessible 
parameter space, the one-loop corrections reach —20% and more for BR(t2 tXk)y ^ = 1)2, 
which can exceed the anticipated ILC precision. In general a strong variation of the one-loop 
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effects with (fAt oii the branching ratios is found, where very large corrections are found in 
S2 for /c = 3,4, where the one-loop contributions can change the BRs by up to —40%. 

The results for T{t2 — )■ bXj) (j = 1, 2) are shown in Figs. [321 ESI For r(t2 — >■ bx^) the 
decay width changes substantially with ^pAf The relative corrections are mostly between 
—5% and —20%, except in S2 for ipAt ~ tt, where r(t2 bxt) becomes very small. The 
absorptive contributions lead to a visible shift in the relative one-loop corrections in S2, 
where the largest effects are found again where T{t2 — )■ hxl) is small. For T{t2 — ?■ ^xt) we 
find a similar size of the corrections. Larger relative corrections of up to —32% are reached 
only in SI where the decay width itself becomes very small. Within SI the larger branching 
ratio values of 7% - 12% are found for t2 — )■ hxt- Here the relative corrections are between 
—18% and —30%, with some variation induced by the absorptive contributions, which can be 
relevant for the LHC and the ILC. For t2 — ?■ hxt larger BR is found in S2, where values 
around 8% are found. The one-loop effects nearly reach —40%, which can be relevant for the 
LHC. Finally it should be noted that the apparently very large corrections on BR(t2 — ^ hxt) 
in S2 (see the lower right plot in Fig. [32]) do not correspond to a negative BR. At ipAt ~ ^ the 
loop corrections are negative and comparably to the (very small) tree-level width, leading 
to BR^"^^ -C BR*'^®^ in Eq. (12081) . The effect of these relatively large loop corrections around 
(pAt ~ can be sizably lowered by including higher-order corrections as, e.g., |A^ioopP- 

Now we turn to the decay modes involving scalar bottom quarks, which have also been 
analyzed in Ref. [38]- In Figs. [341 [35] the results for T{t2 — )■ hH^) {i = 1,2) are presented. 
While we find 1(^2 -5- hH^) at the ~ 1 GeV level, 1(^2 -> hH'^) is only around the 
0.1 GeV level. The relative variation of T(t2 hiH'^) ranges from ~ —13% for large values 
of the width to ~ —23% for small values, with some variations induced by the absorptive 
contributions. For r(i:2 — >■ ^2-^^^) in SI the relative variation can become very large, —60% 
(with a clearly visible shift from the absorptive contributions) but the partial decay width is 
negligibly small. Within S2, where r(t2 h2H^) ~ 0.1 GeV is realized the relative variation 
is at the 10% level. As for the variation with m^-^ we find only small values for the branching 
ratios at the 1% (0.1%) level for the decay to the lighter (heavier) sbottom. The one-loop 
effects on the BRs are only important if other channels are kinematically suppressed. In this 
case the effects can be of the same order as for the partial decay widths itself. Again, the 
apparent very large effect on BR(t2 — > b2H~^) in SI still corresponds to a positive BR; see 
Eq. ([208]). 

The other decay modes involving scalar bottom quarks, t2 — >■ biW^ {i = 1,2) are ana- 
lyzed in Figs. [3S1 [33 As in the analysis with m^^ varied, we find r(t2 — ^ biW~^) at the 
11 (25) GeV level in SI (S2). The relative correction without absorptive contributions 
changes sizably in SI, ranging between 0% and +12%. Taking into account the absorp- 
tive contributions this strongly reduces to 5% and 7%. Within S2 the corrections without 
absorptive contributions are around -|-5% for all ipAt, but the absorptive contributions have 
the opposite effect of strongly enhancing the variation. T{t2 — )■ b2W~^), again as in Sec. 14.21 
is very small and stays below ~ 0.03 GeV. The variation of this negligibly small partial 
decay width is found to be between —6% and +6%; the shift from the absorptive contri- 
butions remains relatively small. Consequently, a relevant branching ratio is found only 
for ^2 biW~^, where values around ~ 20% (~ 16%) are found in SI (S2). The relative 
effects of the one-loop corrections can reach —7% (—10%), which is potentially important 
for physics at the ILC and the LHC. 



56 



r/GeV 




jr/r 


;ree 






10 






22% 








1 













9 

8 


r 


20% 
18% 


\ * 






7 




SI, tree ' 


iD/o 








6 




SI, full . 

SI, abs 


14% 








o 




S2, tree ■ 


12% 






1 


4 




S2 full 

kJ ^ > 1 1111 

S2, abs 


10% 




** 'V ''. 

\ \ 
** '.'V 


/ / 


3 
2 
1 


" 

1 


- 

1 1 1 1 1 1 


8% 
6% 
4% 


_ SI, full 

51, abs 

- S2, full 

52, abs 

1 1 


'\ *• *•. 

'X. \ \ .X, 

V > 

1 1 1 


/ 

/ ,y 



0° 45° 90° 135° 180° 225° 270° 315° 360° 0° 45° 90° 135° 180° 225° 270° 315° 360° 




Figure 23: V{^2 Tree-level ("tree") and full one-loop ("full") corrected partial decay 

widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with LpAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 24: T{i2 — > ^1/12)- Tree-level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 25: T{i2 — > iihs). Tree-level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 26: T{i2 — )■ iiZ). Tree- level ("tree") and full one- loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [H), with ip^t varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 27: T(i2 — tg). Tree-level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [H), with ip^t varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 



61 





Figure 28: T(i2 — > ^Xi)- Tree- level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 29: T(i2 — > ^X^)- Tree- level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 30: T(i2 — > ^Xs)- Tree- level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 31: T(i2 — > txD- Tree-level ("tree") and full one-loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 



65 





Figure 32: T(i2 — )■ bxt)- Tree- level ("tree") and full one- loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 33: T(i2 — > bxt)- Tree- level ("tree") and full one- loop ("full") corrected partial decay 
widths are shown. Also shown are the full one-loop corrected partial decay widths including 
absorptive contributions ("abs"). The parameters are chosen according to SI and S2 (see 
Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the upper right 
plot shows the corresponding relative size of the corrections. The lower left plot shows the 
BR; the lower right plot shows the relative correction of the BR. 
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Figure 34: T(i2 hiH^). Tree- level ("tree") and full one-loop ("full") corrected partial 
decay widths are shown. Also shown are the full one-loop corrected partial decay widths 
including absorptive contributions ("abs"). The parameters are chosen according to SI and 
S2 (see Tab. [I]), with tpAt varied. The upper left plot shows the partial decay width; the 
upper right plot shows the corresponding relative size of the corrections. The lower left plot 
shows the BR; the lower right plot shows the relative correction of the BR. 
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Figure 35: r(t2 — ^ &2-f^^)- Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths are shown. Also shown are the full one-loop corrected partial decay widths 
including absorptive contributions ("abs"). The parameters are chosen according to SI and 
S2 (see Tab. [I]), with LpAt varied. The upper left plot shows the partial decay width; the 
upper right plot shows the corresponding relative size of the corrections. The lower left plot 
shows the BR; the lower right plot shows the relative correction of the BR. 
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Figure 36: T(i2 — )■ biW~^). Tree- level ("tree") and full one-loop ("full") corrected partial 
decay widths are shown. Also shown are the full one-loop corrected partial decay widths 
including absorptive contributions ("abs"). The parameters are chosen according to SI and 
S2 (see Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the 
upper right plot shows the corresponding relative size of the corrections. The lower left plot 
shows the BR; the lower right plot shows the relative correction of the BR. 
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Figure 37: T(i2 — )■ b2W~^). Tree-level ("tree") and full one-loop ("full") corrected partial 
decay widths are shown. Also shown are the full one-loop corrected partial decay widths 
including absorptive contributions ("abs"). The parameters are chosen according to SI and 
S2 (see Tab. [I]), with ipAt varied. The upper left plot shows the partial decay width; the 
upper right plot shows the corresponding relative size of the corrections. The lower left plot 
shows the BR; the lower right plot shows the relative correction of the BR. 
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4.4 The total decay width 

Finally we show the results for the total decay width of t2- In Fig. |38] the upper panels 
show the absolute and relative variation with m^^. The lower panels depict the result for 
varying ^pAf In SI for small m^-^, m^-^ + m^^ ^ 1000 GeV the size of the relative corrections 
of Ftot ranges between +15% and +8%. For larger m^-^ in the two numerical scenarios the 
variation ranges between ~ +7% down to ~ —5% for m^^ = 3 TeV. The variation with ip^t 
is found to be large in both numerical scenarios. Within SI we find values of the relative 
correction between +13% and +7%, decreasing to a range of +9.5% and +11% once the 
absorptive contributions are taken into account. For S2 the absolute values as well as the 
relative correction of Ftot, are larger than in SI. The size of the relative corrections ranges 
between +11% and +15.5%, where the absorptive contributions do not change the overall 
size of the effects but only affect the dependence on tpAt ■ 

5 Conclusions 

We evaluate all partial decay widths corresponding to a two-body decay of the heavy 
scalar top quark in the Minimal Supersymmetric Standard Model with complex parame- 
ters (cMSSM). The decay modes are given in Eqs. ([I]) - ©. The evaluation is based on 
a full one-loop calculation of all decay channels, also including hard QED and QCD radia- 
tion. Such a calculation is necessary to derive a reliable prediction of any two-body decay 
branching ratio. Three-body decay modes can become sizable only if all the two-body de- 
cay channels are kinematically (nearly) closed and have thus been neglected throughout the 
paper. 

We first reviewed the one-loop renormalization procedure of the cMSSM, which is relevant 
for our calculation. This includes the t/t and h/h sector (which has been chosen according to 
the analysis in Ref. |3B]), the gluino sector, and the strong coupling constant. We furthermore 
reviewed the required renormalization of the Higgs and SM gauge boson sector as well as 
the chargino and neutralino sector in the cMSSM. 

We have discussed the calculation of the one-loop diagrams, the treatment of UV and IR 
divergences that are canceled by the inclusion of (hard and soft) QCD and QED radiation. 
Our calculation set up can easily be extended to other two-body decay modes in the cMSSM. 
In fact in order to test our method we checked the finiteness of various other partial decay 
widths (considering neutralino, chargino, and Higgs boson decays). 

For the numerical analysis we have chosen two parameter sets that allow simultaneously 
all two-body decay modes (but could potentially be in conflict with the most recent SUSY 
search results from the LHC). The masses of the scalar top quarks in these scenarios are 
260 and 650 GeV, and 720 and 1200 GeV for the lighter and the heavier stop, respectively. 
Consequently, both scenarios result in copious scalar top quark production at the LHC. A 
decay of the heavy stop to a lighter stop (or sbottom) and a neutral (or charged) Higgs 
boson can serve as a source of Higgs bosons at the LHC, thus a precise knowledge of stop 
branching ratios is required. The first scenario also allows i\i2 production at the ILC(IOOO), 
where statistically dominated experimental measurements of the heavy stop branching ratios 
will be possible. Depending on the integrated luminosity a precision at the few percent level 
seems to be achievable. 
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Fi gure 38: Fjot- The tree-level ("tree") and full one-loop ("full") corrected total decay widths 
shown with the parameters chosen according to SI and S2 (see Tab. [1]). The upper left plot 
shows the total decay width; the upper right plot shows the corresponding relative size of the 
total corrections, with m^-^ varied. The lower plots show the same but with ipAt varied. Also 
shown is the full one-loop corrected total decay width including absorptive contributions 
("abs"). 
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In our numerical analysis we have shown results for varying m^^ and (fAt, the phase 
of the trilinear coupling At. In the results with varied m^^ only the lighter values allow 
i\i2 production at the ILC(IOOO), whereas the results with varied ipAt have sufficiently light 
scalar top quarks to permit e~^e~ — >■ i\i2- In the two numerical scenarios we compared the 
tree-level partial widths with the one-loop corrected partial decay widths. In the analysis 
with (fiAt varied we showed explicitly the effect of the absorptive parts of self-energy type 
corrections on external legs. We also analyzed the relative change of the partial decay widths 
to demonstrate the size of the loop corrections on each individual channel. In order to see the 
effect on the experimentally accessible quantities we also show the various branching ratios 
at tree-level (all channels are evaluated at tree-level) and at the one-loop level (with all 
channels evaluated including the full one-loop contributions). Furthermore we presented the 
relative change of the BRs that can directly be compared with the anticipated experimental 
accuracy. 

We found sizable, roughly (9(10%), corrections in all the channels. For some parts of the 
parameter space (not only close to thresholds) also larger corrections up to 30% or 40% have 
been observed. This applies especially to the BR(t2 tihn) with n = 1, 2, 3. The size of the 
full one-loop corrections to the partial decay widths and the branching ratios also depends 
strongly on ipAf The one-loop contributions, again being roughly of (9(10%), often vary by 
a factor of 2 — 3 as a function of ipAt ■ In some cases the absorptive contributions can change 
the result visibly. All results are given in detail in Sees. 14.21 and 14.31 

The numerical results we have shown are, of course, dependent on the choice of the 
SUSY parameters. Nevertheless, they give an idea of the relevance of the full one-loop 
corrections. The largest partial decay width, if kinematically allowed, is T{t2 — )■ tg) in our 
scenarios, dominating the total decay width, Ftot, and thus the various branching ratios. For 
other choices of with rrig > mi^ the corrections to the partial decay widths would stay 
the same, but the branching ratios would look very different. Decay channels (and their 
respective one-loop corrections) that may look unobservable due to the smallness of their 
BR in our numerical examples could become important if other channels are kinematically 
forbidden. 

Following our analysis it is evident that the full one-loop corrections are mandatory for 
a precise prediction of the various branching ratios. The results for the scalar top quark 
decays will be implemented into the Fortran code FeynHiggs. 
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Appendix: Absorptive parts from self-energy type con- 
tributions 



As indicated in the main text, contributions to the partial decay widths can arise from the 
product of the imaginary parts of the loop functions (absorptive parts) of the self-energy 
type contributions in the external legs and the imaginary parts of complex couplings enter- 
ing the decay vertex or the self-energies. In our calculation these corrections are taken into 
account via wave function correction factors bZ (which should not be confused with the field 
renormalization constants bZ^ which have been introduced via the multiplicative renormal- 
ization procedure). For the off-diagonal wave function correction factors, this procedure has 
been checked against explicitly including the (renormalized) self-energy type corrections of 
the external legs, and full agreement has been found. The corrections from the absorptive 
parts can be sizable. 

It is possible to combine the wave function correction factors with the field renormaliza- 
tion constants in a single Z factor, Z\ see e.g. Ref. [92] and references therein. 

However, if the external particles were stable the wave function corrections would be fully 
taken into account via the field renormalization constants. The Z factors listed in Sec. [2] 
also ensure that the external (stable) particle does not mix with other fields, which is one of 
the on-shell properties. In our scenarios, this is true e.g. the lightest neutralino. 

In the case of quasi stable particleqfj, additional contributions to the mixing can occur 
so that the field renormalization constants only partly ensure no mixing (for a more detailed 
explanation; see the subsection scalar quarks below). Extra diagonal contributions can also 
be taken into account via Z factorsjf^ Here we briefly list all the resulting constants. 

Scalar quarks 

(i) For an on-shell particle state, no mixing with another state should occur, corresponding 
to 

S^.,(m|) = , S,-,,(m|j = , S,-,,(m|) = , t-,,,[m\^ = , (214) 

in the case of squarks. Partly, this is already fulfilled due to our renormalization 
conditions Eq. (159|) in Sec. 12.1.11 As the considered external particles are quasi stable, 
in spite of the renormalization conditions above, there remains a contribution of the 
imaginary parts of the loop functions. This contribution is taken into account via 
wave function correction factors 6Z, which are different for incoming squarks/out going 
antisquarks (unbarred) and outgoing squarks/incoming antisquarks (barred), 

2^^^i?i^ , (215) 
{ml - ml) 

^'^ Which means that considering them as external particles is an approximation, which is justified because 
in our decays the contributions from the (additional) diagonal SZ are numerically rather negligible. 

There is still an ongoing discussion whether the diagonal field renormalization constants take into 
account all the contributions needed to ensure the on-shell properties of the external particles or whether an 
extra wave function correction factor 6Z is needed. 
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where Im takes only the imaginary part of the loop functions. Compact expressions 
for practical numerical calculations are obtained via the combined Z factors 5Z, 
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(ii) The diagonal contributions result in the following combined Z factors: 



Res'; (m? 



zImS'r.fm? 



[^^d.= [^Z, + 5Z,],^^ = 



(217) 
(218) 
(219) 
(220) 

(221) 
(222) 



It should be noted that ReSij(/) = (ReSji(p2))* holds due to CVT invariance and the 
squark field renormalization constants obey [Re5Zg]jj = [Re5Zg]*j. = [5Zq]*j. , which is 
exactly the case without absorptive contributions as described in Sec. 12.1.11 

In the following we will only give the Z factors that combine the renormalization factors 
and the additional wave function correction factors. The derivation is analogous to the one 
performed in the squark sector. 



Quarks 

The new (diagonal) combined factors Z^^ taking into account the absorptive part of the self- 
energy type contribution on an external quark leg are different for incoming quarks/outgoing 
antiquarks (unbarred) and outgoing quarks/incoming antiquarks (barred), 



^'""iO + ml { E^'iml) + Sf (m^)) + ( Ef' (m^) + 



± 



2mr. 



fK) 



Sj/^(m5)+m^fs^'(m2^ 



ml) + (mj; 



2 m, 



f(mJ)-Sf(mJ 



(223) 



(224) 



The diagonal quark field renormalization constants obey KeSZq^^ = [dZq^^]*, which is 
exactly the case without absorptive contributions as described in Sec. 12.1.11 

There are no additional off-diagonal terms to the absorptive contributions because the 
CKM matrix has been set to unity. 
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Gluinos 

The new combined factors Zg^ taking into account the absorptive part of the self-energy type 
contribution on the external gluino leg are unbarred (barred) for an incoming (outgoing) 
gluino, 




The last formula holds due to the Majorana character of the gluino and the Z~g factors obey 
Re 5Zg = Re SZg = SZ^ SZg, which is exactly the case without absorptive contributions 
as described in Sec. I2.1.2I 



Higgs bosons 

Finite contributions from the neutral Higgs wave function correction factors are taken into 
account via the Z matrix; see Eq. fll32p . which is a complex quantity. The application of 
the Z matrix at the amplitude level automatically takes any absorptive contribution into 
account. 

For the charged Higgs bosons, the new combined factors Zh-h+ (unbarred (barred) for 
an incoming (outgoing) Higgs) read 

6Zh-h+ = -^'H-H+iMl±) , (227) 
5Zh-h+ = 5Zh-h+ (228) 



instead of Eq. ( |99gP in addition with Eq. f ll34p . 



Vector bosons 

For the vector bosons, the new combined factors Zs^yY^rzz} are 

5Zzz = -^I'ziMl) , 6Zww = -^^wiM^) , (229) 

SZzz = SZzz , SZww = S2ww ■ (230) 

However, we found that the additional corrections from vector boson self-energies due to the 
imaginary parts do not give a contribution (because in this paper all SUSY masses are larger 
than Mz), and hence no change in the Z factors is required. 



Charginos and neutralinos 

More details to the new combined factors Z^- and Z^o (taking into account the absorptive 
part of the self-energy type contributions on the external legs) can be found in Ref . [53] • In 
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our notation they read (unbarred for an incoming neutralino or a negative chargino, barred 
for an outgoing neutralino or negative chargino) 
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(237) 
(238) 



The chargino/neutrahno Z factors obey Re SZ^^^ = [Re dZ^^^^ = [6Z'^^^]\ which is exactly 

the case without absorptive contributions as described in Sec. 12.31 or in other words ^Z^^^ = 

[5Zj/^]t + [5Z|/^]. The Eqs. m7\f and mE\f hold due to the Majorana character of the 
neutralinos. 
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